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The accelerated expansion of space during the period of cosmological inflation leads to trans- 
Planckian issues which need to be addressed. Most importantly, the physical wavelength of fluctua- 
tions which are studied at the present time by means of cosmological observations may well originate 
with a wavelength smaller than the Planck length at the beginning of the inflationary phase. Thus, 
questions arise as to whether the usual predictions of inflationary cosmology are robust considering 
our ignorance of physics on trans-Planckian scales, and whether the imprints of Planck-scale physics 
are at the present time observable. These and other related questions are reviewed in this article. 
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^ . I. INTRODUCTION 

The theory of cosmological inflation [1] (see als [2-4] for related original work) has become the current paradigm 
of early universe cosmology. Not only does the scenario solve some of the conceptual problems of Standard Big Bang 
cosmology, the previous paradigm of cosmology, but it provided [5] the first explanation for the origin of the large-scale 
structure of the universe based on causal physics (see also [2, 6, 7] for related original work). Better yet, inflationary 
cosmology was predictive. The scenario predicted an almost scale-invariant spectrum of curvature fluctuations which 
are coherent and passive on scales which in the early universe are larger than the Hubble radius. As has been known 
since around 1970 [8, 9], such a spectrum predicts characteristic oscillations in the angular power spectrum of cosmic 
microwave (CMB) anisotropy maps, oscillations which were first discovered by the Boomerang experiment [10] and 
later confirmed with high accuracy by the WMAP satellite experiment [11]. 

Inflationary cosmology is based on the assumption that there is a period in the very early universe during which 
space is expanding at an accelerated rate (typically almost exponentially). Figure 1 depicts a space-time sketch of 
inflationary cosmology. The vertical axis is time, the horizontal axis represents physical distance. The period of 
inflation begins at a time t\ and ends at time t^. The accelerated expansion of space during this time interval can 
create a large and almost spatially flat universe from an initial state in which the size of space is microscopic. In 
this way, inflation solves the "size problem" of Standard Cosmology. In an accelerating universe the contribution of 
spatial curvature to the total density decreases. Hence, a spatially flat universe is a local attractor in initial condition 
space [12] l , thus providing a solution to the "flatness problem". 

In Figure 1 three distance scales are shown. The solid curve which corresponds to an almost constant value during 
the period of inflation is the Hubble radius £n(t) which is given by the inverse expansion rate of space. The Hubble 
radius plays an important role in the evolution of cosmological perturbations, the inhomogeneities which give rise 
to the large-scale structure of the universe and which lead to CMB anisotropics . On length scales smaller than the 
Hubble radius, matter fluctuations oscillate almost like in flat space-time. These fluctuations freeze out when the 
wavelength crosses the Hubble radius. On larger scales, the evolution of the perturbations is governed by gravity. 
The dashed curve which tracks the Hubble radius before the onset of inflation but whose length grows exponentially 
during the period of inflation is the horizon, the forward light cone from a point at the time of the Big Bang. The 
horizon is the largest distance that causal information can travel. The fact that the horizon becomes exponentially 
larger than the Hubble radius provides a solution to the "horizon (or homogeneity) problem" of Standard Cosmology: 
In the context of inflation, the causal horizon can be much larger that any scale which we observe today. In particular 
this can explain the near isotropy of the CMB. 
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1 Note that it is not a global attractor, and hence a certain degree of spatial flatness is required in order to be able to enter into a period 
of accelerated expansion. 
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FIG. 1: Space-time sketch of inflationary cosmology. The vertical axis is time, the horizontal axis corresponds to physical 
distance. The solid line labelled k is the physical length of a fixed comoving fluctuation scale. The role of the Hubble radius 
and the horizon are discussed in the text. 

The third length scale shown in Figure 1, the curve labelled by k, represents the physical length of a fluctuation 
mode. The fact that modes well inside the Hubble radius at early times become exponentially larger than the Hubble 
radius by the end of the period of inflation provides a causal mechanism by which fluctuations are generated in 
inflationary cosmology. Since classical matter redshifts exponentially during the period of inflation, it is natural 
to assume that matter approaches the vacuum state. The vacuum, however, is permeated by quantum vacuum 
perturbations, and hence it was conjectured [2, 6, 7] that the primordial fluctuations are quantum mechanical in 
nature. The theory of structure formation in inflationary cosmology is based on the quantum theory of cosmological 
perturbations [13, 14] (see e.g. [15, 16] for in-depth reviews, and [17] for a pedagogical overview). Based on the above 
discussion it is reasonable to assume that fluctuations begin in their quantum vacuum state [5] . While on sub- Hubble 
scales, the fluctuation modes undergo quantum vacuum oscillations. They freeze out once the wavelength of the 
mode becomes larger than the Hubble radius. Thereafter the modes are squeezed while on super-Hubble scales. In 
particular, all modes acquire the same angle in phase space. They hence enter the Hubble radius at late times as 
standing waves, leading to the characteristic acoustic oscillations in the angular power spectrum of CMB anisotropies 
which were mentioned at the beginning of this section. Based on the time-translation symmetry of the inflationary 
phase, it is expected [6] that the spectrum of the density perturbations produced during inflation is approximately 
scale-invariant, and the mathematical analysis confirms this result. 

The causal generation mechanism of fluctuations described above is the most significant success of inflationary 
cosmology. However, a second look at the space-time sketch of Fig. 1 reveals a serious conceptual problem [18, 19]. 
A sufficiently long period of inflation is required in order for modes which are currently probed in cosmological 
observations to have a wavelength smaller than the Hubble radius at the beginning of the period of inflation - which 
is a necessary condition for the validity of the inflationary structure formation scenario. However, if the period of 
inflation is only slightly longer (70 e-foldings in models in which the energy scale at which inflation takes place is 
close to the scale of Grand Unification), then the wavelengths of all fluctuation modes which are currently inside the 
Hubble radius were smaller than the Planck length at the beginning of the period of inflation. We do not understand 
physics on length scales smaller than the Planck length. The problem for inflationary cosmology is that fluctuation 
modes emerge from this sub-Planck-wavelength zone of ignorance, as is sketched in Fig. 2. New physics is required 
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FIG. 2: Space-time diagram (sketch) of inflationary cosmology where we have added an extra length scale, namely the Planck 
length £ pl (majenta vertical line). The symbols have the same meaning as in Figure 1. Note, specifically, that - as long as the 
period of inflation lasts a couple of e-foldings longer than the minimal value required for inflation to address the problems of 
standard big bang cosmology - all wavelengths of cosmological interest to us today start out at the beginning of the period of 
inflation with a wavelength which is smaller than the Planck length. 

to truly understand the origin and early evolution of the fluctuations 2 . 

Based on the geometry of Figure 2, there are good reasons to expect that Planck-scale physics can lead to a 
modification of the spectrum of cosmological perturbations: If wc imagine starting off all perturbation modes on a 
fixed space-like Cauchy surface at the beginning of the period of inflation, then short distance modes will feel the 
effects of the modified physics for a longer time than long wavelength modes. If the evolution of the modes is not 
adiabatic on short wavelength scales, then one would expect the spectrum for short wavelength modes to be boosted 
relative to that of long wavelength modes. This would produce a blue tilt in the spectrum of fluctuations. 

The above problem is now called the trans- Planckian problem for inflationary cosmology. While it may be a problem 
from the point of view of inflationary theory, it can instead to viewed as a trans- Planckian window of opportunity to 
probe Planck-scale physics with current cosmological observations. The accelerated expansion of space brings scales 
into the observable range which we have no chance of probing in accelerator physics. 

In the following sections, we will review these trans-Planckian issues for inflationary cosmology. In Section 2 we 
discuss some theoretical approaches, beginning with a brief review of some of the original approaches to studying the 
sensitivity of the predictions of inflationary cosmology to assumptions about trans-Planckian physics. We continue 
with a discussion of two more recent realizations of inflation in the context of specific models of trans-Planckian 
physics. In one case we find dramatic differences compared to the "standard results" , in the second case we recover 
the usual spectrum. We end Section 2 with some general remarks. In Sections 3 and 4 we discuss observational 
constraints on the magnitude of trans-Planckian effects in inflationary cosmology. Finally, in Section 5 we attempt to 
draw connections with some broader issues. 



2 As stressed in [20] in a different context, for a model of structure formation to be consistent, it is not sufficient that the energy scale 
of the background cosmology is smaller than the Planck scale. If the wavelength of the fluctuation modes is smaller than the Planck 
length, the predictions of the model cannot be trusted. 
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II. THEORETICAL APPROACHES 



A. Original Approaches 

1. Modified Dispersion Relations 

The original work on the trans-Planckian problem for inflation was performed in the context of modified dispersion 
relations, [19, 21] for the linear fluctuation modes, following works [22, 23] which studied the dependence of black 
hole radiation on Planck scale physics. This method represents a toy model to study the unkown effects of trans- 
Planckian physics, a method which has the practical advantage of remaining within the context of applicability of 
linear cosmological perturbation theory. Before introducing the basics of the method, we must remind the reader of 
the relevant formalism. 

We are interested in following the scalar metric perturbations, the linearized gravitational fluctuation modes induced 
by matter perturbations. In the context of General Relativity, The relevant canonical variable (the Mukhanov-Sasaki 
variable [13, 14]) v satisfies the following Fourier space equation of motion [15, 17] 
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where k denotes the comoving momentum and a prime stands for the derivative with respect to conformal time 77. The 
function z(rj) depends on the background cosmology. If the equation of state of the background is time- independent, 
then z(r]) is proportional to the scale factor 0(77). The variable v is related to the curvature fluctuation £ in comoving 
gauge, a coordinate system in which the matter fluctuation vanishes, via £ = z~ 1 v. The variable £ carries vanishing 
mass dimension. 

Of particular interest is the dimensionless power spectrum V^(k) defined via [for a more accurate definition, see 
Eq. (23)] 



(2) 



The power spectrum is called "scale-invariant" if V^(k) is independent of k. 

The idea of the modified dispersion relation method is to introduce a non-trivial relation between the physical 
frequency and momentum of the modes 
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where the function L} p h ys (u) deviates from its usual form w p hy S ('u) = u only in the ultraviolet (UV). As will be discussed 
below, such dispersion relations arise from specific theories of Planck-scale physics such as Hofava-Lifshitz gravity 
[24]. 

The modification of the dispersion relation implies a change in the Hamiltonian of the perturbation modes which 
has important effects at short wavelengths. The analysis of [19] is based on the following assumptions: first, we 
consider a space-like hypersurface on which we set up initial conditions. Second, the initial conditions arc chosen to 
represent the lowest energy state of the local Hamiltonian at the initial time. Since this Hamiltonian will deviate 
from the usual one on short wavelength scales, the initial state will deviate from the usual Bunch-Davics vacuum of 
the un- modified theory, the deviations increasing towards the ultraviolet. 

For mild distortions of the dispersion relation such as those considered in [22] where u> 2 asymptotes to a constant 
in the ultraviolet, the evolution of the mode functions will be adiabatic (they are the WKB solutions) and track the 
local vacuum state. In this case, the evolution will lead to the usual scale-invariant spectrum at late times. 

On the other hand, if the dispersion relation violates the adiabaticity condition 
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then an initial vacuum state will evolve to a state which on large scales differs from the Bunch-Davies vacuum of 
the unmodified theory. The deviations increase as the frequency increases, thus leading to a steep blue spectrum in 
violation of the observational constraints (see Sections 3 and 4). 

Deviations from adiabaticity arise for one branch of the modified dispersion relations proposed in [23] . Specifically, 
one can consider the example [25, 26] 
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with bi and 62 being two positive constants. In the case of this dispersion relation (sketched in Fig. 7) there can be a 
region of values of £; p hy S for which the adiabaticity condition is violated. In this example, adiabaticity is maintained 
both in the far ultraviolet and the infrared. 

While we do not expect concrete models of modified physics on trans-Planckian scales to yield dispersion relations 
which are adiabatic in the far UV, assuming that this takes place in our toy model makes it easy to justify our 
initial conditions: we can start modes in the adiabatic vacuum in the far UV. As the universe expands, the physical 
wavenumber will redshift. Modes which begin in the far UV adiabatic regime will experience a time interval of finite 
duration during which the wavenumber is in the non-adiabatic region. Hence, the wavefunction will obtain corrections 
from its WKB form. In particular, the amplitude of the final fluctuations will change [25]. 

To a first approximation (in the deviation of the expansion of space from exponential), the shape of the spectrum 
will not change since all modes spend the same amount of time in the region of non-adiabaticity (as discussed in 
Section 3, there will be superimposed oscillations in the spectrum). However, in models in which non-adiabaticty is 
violated at all UV scales (like the ones considered in [19]), short wavelength modes spend more time in the region of 
non-adiabaticity, and hence the spectrum of cosmological perturbations will acquire a blue tilt whose spectral slope 
can well exceed current limits. 



2. New Physics Hypersurface 

A more conservative approach to the trans-Planckian problem is simply not to evolve the fluctuation modes during 
the time period in which their wavelength is smaller than the length scale of the new physics. This corresponds to 
introducing a time-like new physics hypersurface on which initial conditions arc imposed. We expect the wavelength 
at which the new physics hypersurface is reached to be given by the Planck length (or the string length - in the 
context of string cosmology - if it is larger) . 

At this point, the trans-Planckian problem has simply been shifted to the problem of choosing initial conditions on 
the new physics hypersurface. The most conservative approach is to start modes off in their local adiabatic vacuum. 
Note, however, that precisely the main point of the analysis of the previous subsection is that new physics acting on 
trans-Planckian scales can well produce states which are very different from the local adiabatic vacuum at the time 
when the wavelength equals the new physics scale. 

Even continuing with the conservative approach mentioned above, a careful tracking of the mode wavefunctions 
reveals residual oscillations in the power spectrum [27-31]. The amplitude of these oscillations depends on whether 
we are computing the spectrum of test scalar fields on a fixed background metric, gravitational waves, or scalar metric 
fluctuations. In the two former cases the amplitude of the oscillations in the power spectrum is 

H ,2 



A GW ~ 0(1)1 _j , (6) 

where M c is the mass scale of the new physics hypersurface, and in the case of scalar cosmological perturbations we 
have 

*~0(i>£. CD 

Note that within the new physics hypersurface approach, the setting of the initial state is not unique. Whereas the 
prescription adopted by [27, 28, 31] leads to oscillations of the abovementioned magnitude, the choice made by [29] 
and studied in more detail in [32] leads to oscillations with amplitude suppressed by three powers of H/M c . The 
dependence of the amplitude of the trans-Planckian effects on the definition of the "local vacuum state" was explored 
in detail in [30] . 

In the context of the choice of the vacuum state for cosmological perturbations there has been some discussion about 
vacuum states which are alternatives to the usual Bunch-Davies [33] vacuum. These are the a vacua [34] discussed 
in the context of cosmological perturbations in [28, 35]. They are in different superselection sectors of states than the 
Bunch-Davies vacuum, and applied to inflationary cosmology they also lead to oscillations in the power spectrum of 
fluctuations. The use of such vacua has, however, been criticized on the basis of their instabilities [36, 37] (but see 
[38]). Our view, however, is that in the context of inflationary cosmology which is not past eternal [39], it makes more 
sense to restrict attention to states which are in the same superselection sector as the usual vacuum, and which can 
hence be generated from such a vacuum by local trans-Planckian physics. 
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3. Other Approaches 

In the two previous subsections we discussed two general frameworks for studying trans-Planckian effects on the 
spectrum of cosmological perturbations. Another general framework is the effective field theory approach. The idea 
of effective field theory is to integrate out the high energy physics and to derive an effective Lagrangian which only 
involves observable scales. The application of effective field theory techniques to the trans-Planckian problem of 
inflation was pioneered in [40] and [41]. However, the expansion of space (in particular the fact that Planck-scale 
modes are redshifted into the visible sector) leads to challenges which are not present in usual applications of effective 
field theory 3 , challenges which were not taken into account in [40], where it was claimed that trans-Planckian effects 
on the spectrum of cosmological perturbations cannot be larger than 0(H/M c ) 2 . However, as shown e.g. in the work 
of [41], even in the context of effective field theory and local Lorentz invariance it is possible for the pre-inflationary 
dynamics to produce states which deviate from the standard Bunch-Davics vacuum and for which hence the corrections 
to the power spectrum are larger. 

Some authors have considered general modifications of physics on Planck scales such as modified uncertainty prin- 
ciples [43], space-space non-commutativty [27, 44] and space-time non-commutativity [45] and computed the changes 
in the spectrum of perturbations in an inflationary cosmology. All of these approaches can be viewed as leading to 
specific prescriptions of setting the intial conditions on a space-like new physics hypcrsurface. 

If we make the strong assumptions firstly that the microscopic structure of space-time on Planck scales is Lorentz 
invariant and secondly that the state of the system is in its local vacuum, then it indeed follows that the corrections 
to the spatial correlation function and hence to the power spectrum are bounded by 0(H/M c ) 2 [46]. Even if one 
accepts the first assumption, the second one is a very strong one. We know that (in the context of Einstein gravity 
with matter fields obeying the usual energy conditions) inflation is singular in the past [39], and there is thus pre- 
inflationary dynamics which must be taken into account when determining what the initial state of the fluctuation 
modes is. In particular, perfectly Lorentz-invariant theories can lead to a bouncing universe, e.g. through the addition 
of quintom matter [47] or in the context of specific string theory backgrounds [48] . 

B. Specific Models 

1. Bouncing Inflation 

In theories which admit non-singular bouncing cosmologies (e.g. Hofava-Lifshitz gravity in the presence of non- 
vanishing spatial curvature [49]), it is possible to have a post-bounce inflationary phase and, by time reflection 
symmetry, a contracting deflationary phase (see Figure 3 for the corresponding space-time diagram). The spectrum 
of cosmological perturbations was studied in this context in [50] . 

Since the universe begins large and cold, it is natural to start with perturbations in their Bunch-Davies vacuum 
state. In a contracting universe the curvature fluctuation variable £ grows on super-Hubble scales (see e.g. [51, 52]). 
In a radiation phase of contraction the growth rate is proportional to r\~ x (the conformal time i] is increasing towards 
0), and in the deflationary phase £ increases as rf (n is increasing to oo in this phase). In this process, the initial 
vacuum spectrum is transformed into a dimcnsionless spectrum which scales as 

V[k,-U(k)} ~ k~\ (8) 

where — ij(fc) is the time when the scale k enters the Hubble radius during the deflationary phase. 

The evolution of the fluctuations between the time U(k) of Hubble radius entry during the contracting phase and 
Hubble radius exit at the time U(k) is symmetric. Hence, the spectrum at Hubble radius exit is the same as at Hubble 
radius entrance. Since £ is conserved on super-Hubble scales in the expanding phase, the final spectrum of curvature 
fluctuations in the expanding phase is given by (8), i.e. it is an n s = —3 spectrum. This is a highly red spectrum 
compared to a scale- invariant one (n s = 1 which is indicated by experiments) . The redness of the spectrum is due to 
the extra growth which long wavelength modes undergo since they spend a longer time on super-Hubble scales. 

As will be discussed in the following section, a red spectrum is not only constrained by observations, but also by 
back-reaction. Applied to our model, this implies that the bouncing model cannot be symmetric about the bounce 
point. The phase of deflation must be short. As a consequence, one obtains a scale-invariant spectrum in the ultraviolet 
and a red spectrum in the infrared. 



These challenges also arise in some non-gravitational contexts, e.g. when studying quantum field theory in strong external electromagnetic 
fields [42]. 
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FIG. 3: Space-time sketch of the bouncing inflation model. The vertical axis is time, the horizontal axis denotes the physical 
space coordinate. The Hubble radius H -1 is constant during the deflationary phase of contraction and during the inflationary 
phase of expansion. It diverges around the bounce point. The solid (blue) curve labelled k indicates the wavelength of a 
fluctuation mode which exits the Hubble radius in its Bunch-Davies vacuum state during the radiation phase of contraction 
at time ti(fc), evolves on super-Hubble scales until it enters the Hubble radius at time t3(k) = —ti(k) during the deflationary 
phase. The time ti is the transition between radiation contraction and deflation. 

This model is an example of how trans-Planckian effects can lead to very large deviations from scale-invariance. 
As in most non-singular bouncing cosmologies, it requires specially tuned matter and initial condition modeling to 
obtain a scenario in which a period of deflation follows after a radiation phase. One way, as explored in the quintom 
matter bounce scenarios [47] is to have the pre-bounce radiation come from a scalar field (p with quartic potential. 
While the scalar field is oscillating about its ground state, the time-averaged equation of state is that of radiation. 
In the contracting phase the amplitude of oscillations increases. But once the amplitude exceeds the Planck scale, a 
deflationary slow-climb phase sets in. 

2. Hofava-Lifshitz Inflation 

Hofava-Lifshitz (HL) gravity is a four space-time-dimensional theory which has been proposed [24] as a power- 
counting renormalizable theory of quantum gravity. In this theory, the microscopic Lagrangian is not locally Lorentz- 
invariant. There is a distinguished time direction, and the gravitational Lagrangian contains higher space-derivative 
terms. As a consequence, space-time diffcomorphism invariancc is also lost. The theory is only invariant under the 
subgroup of spatial diffcomorphisms, and under separate space-independent time reparamctrizations. The guiding 
principle of the construction of the Lagrangian is invariance with respect to the residual symmetries and power- 
counting renormalizability with respect to an anisotropic scaling 

x -> bx , t^b 3 t, (9) 

where b is some constant. There are two versions of the theory - the projectable version in which the lapse function 
N 4 is a function of only time, and the non-projectable version in which N can depend on space and time. 



4 In the Hamiltonian approach to gravity, the lapse function is the fluctuation in the time-time component of the metric. 
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Renormalizability allows up to six-space-derivative terms in the equations of motion. Hence, quite naturally a mod- 
ified dispersion relation for fluctuation modes results. There is a second major difference in the theory of cosmological 
perturbations in HL gravity compared to Einstein gravity: because of the reduced symmetry, there is an extra scalar 
gravitational fluctuation mode (see e.g. [53] for a review of HL gravity). It has been shown [54] that in the case 
of the projectable version of HL gravity the extra scalar mode is sick (either ghost-like or tachyonic), while in the 
non-projectable version it can be [55] well-behaved (for suitable choices of the free parameters of the Lagrangian). 
The extra mode in fact becomes massive in the infrared and thus only plays an important role in the ultraviolet. 

Because of the existence of the extra degree of freedom and because of the modified dispersion relation one might 
expect that the spectrum of cosmological perturbations in inflationary cosmology in the context of HL gravity would 
lead to a very different spectrum of cosmological perturbations than in Einstein gravity. However, a careful analysis 
of this issue [56] has shown firstly that the evolution of the regular fluctuation mode is adiabatic and secondly that the 
fluctuations induced by the extra degree of freedom are sub-dominant. Hence, up to oscillations in the spectrum whose 
frequency is too high to be observationally detectable, the spectrum of cosmological perturbations which emerges when 
starting all modes off at a fixed time in their instantaneous minimum energy state is scale invariant. 

This example indicates a certain degree of robustness of the standard predictions of inflation even if the underlying 
physics at the Planck scale is substantially modified. 

C. General Comments 

It is interesting to ask to which extend the trans-Planckian problem discussed here is specific to inflationary 
cosmology or what aspects of the problem arise in any cosmological background. To address this question, let us 
introduce two paradigms alternative to inflation which can explain the observed almost scale-invariant spectrum of 
cosmological fluctuations. The first is the string gas cosmology model introduced as a background cosmology in [57] 
(see also [58]) and shown to lead to a mechanism for producing an almost scale- invariant spectrum of cosmological 
perturbations in [59, 60]. The second is the matter bounce scenario of [51, 52], a non-singular bouncing cosmology 
with an initial matter-dominated phase of contraction. 

String gas cosmology (see e.g. [61] for a recent review) is based on the realization that a gas of closed strings 
has a maximal temperature, the "Hagedorn temperature" [62] Tr. At temperatures close to Xh, all of the modes 
of strings including winding modes are excited. This high temperature phase is called the "Hagedorn phase". It 
is then reasonable to assume that the universe begins in the Hagedorn phase. Since the temperature of a box of 
strings is independent of the radius R of the box for a wide range of values of R (assuming that the entropy of the 
system is large), it is not unreasonable to assume that the Hagedorn phase is quasi-static (including fixed dilaton). 
Eventually (through the decay of string winding modes into string loops) the universe transits into the radiation 
phase of Standard Cosmology. The time when this happens is called £r in analogy to the reheating time at the end 
of inflation. In the context of this string gas cosmology background, it was realized in [59, 60] that thermal string 
fluctuations in the Hagedorn phase induce a scale-invariant spectrum of curvature fluctuations at late times with a 
slight red tilt, and a scale- invariant spectrum of gravitational waves with a slight blue tilt [63] - a key prediction of 
string gas cosmology with which it can be observationally distinguished from inflationary cosmology. 

Figure 4 shows a space-time sketch of string gas cosmology. If the Hagedorn temperature is taken to be similar 
to the post-inflation temperature (which is typically comparable to the scale of Grand Unification), then the post- 
reheating evolution in inflationary cosmology is identical to the post-Hagedorn dynamics in string gas cosmology. 
Working backwards from the present time, it is not hard to see that the physical wavelengths of fluctuations which 
are observed today are many orders of magnitude larger than the Planck scale at the time t^ (1 mm is a typical 
value) . 

Since in string gas cosmology scales which are observed now never had a wavelength close to the Planck scale, the 
trans-Planckian problem for cosmological perturbations is absent. The basic problems of quantum field theory in an 
expanding universe (e.g. time-dependence of the Hilbert space of Fourier modes [64]) still are present, but they are 
not directly coupled to observations, unlike in inflationary cosmology where they are. 

Figure 5 represents a space-time sketch of a matter bounce cosmology. The simplest way to visualize the background 
space-time is to take a mirror inverse of our expanding Standard Cosmology space-time (which is a contracting universe 
which begins in a phase of matter-domination) and match it to the expanding phases of Standard Cosmology via a 
short non-singular bounce phase (see [65] for a recent review of the matter bounce scenario). It goes without saying 
that new physics (either in the matter or in the gravitational sector) is required to obtain such a non-singular bounce. 
If the energy scale of the non-singular bounce is comparable to the scale of Grand Unification, then - as is the case 
in string gas cosmology - the physical wavelengths of scales which are currently probed in cosmological observations 
never enter the sub-Planck length zone of ignorance, and hence there is no trans-Planckian problem for cosmological 
perturbations. 
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FIG. 4: Space-time diagram (sketch) showing the evolution of fixed co-moving scales in string gas cosmology. The vertical axis 
is time, the horizontal axis is physical distance. The solid curve represents the Einstein frame Hubble radius H^ 1 which shrinks 
abruptly to a micro-physical scale at irt and then increases linearly in time for t > t^,. Fixed co- moving scales (the dotted lines 
labeled by fci and /C2) which are currently probed in cosmological observations have wavelengths which are smaller than the 
Hubble radius before £r. They exit the Hubble radius at times ti(fc) just prior to 4r, and propagate with a wavelength larger 
than the Hubble radius until they reenter the Hubble radius at times tf(fc). 
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FIG. 5: Space-time sketch in the matter bounce scenario. The vertical axis is conformal time r\, the horizontal axis denotes a 
co-moving space coordinate. Also, T-L~ x denotes the co-moving Hubble radius. 
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As first realized in [51, 52], initial vacuum fluctuations which exit the Hubble radius during the matter-dominated 
phase of contraction acquire a scale-invariant spectrum of curvature perturbations on supcr-Hubble scales. The key 
point is that curvature perturbations in a contracting universe grow on super-Hubble scales, and in a matter-dominated 
background the extra growth which long wavelength modes experience because they are super-Hubble for a longer 
time is exactly the right amount to convert a blue vacuum spectrum into a scale-invariant one. One of the key 
predictions of the matter bounce scenario is that there is a cutoff scale set by the transition time between matter 
domination and radiation domination during the contracting period below which the spectrum of perturbations turns 
blue [66] . In addition, the growth of the curvature fluctuation on super-Hubble scales leads to a particular shape and 
reasonably large amplitude of the bispectrum [67]. 



III. THEORETICAL FRAMEWORK 



We have seen in the last sections that trans-Planckian physics could modify the inflationary power spectrum. Since 
we measure the inflationary power spectrum when we measure the CMB anisotropics, this opens the possibility to test 
these effects with astrophysical data. In order to carry out this program, we must first understand how to calculate 
the trans-Planckian corrections. At this point, we meet a first difficulty. Since the physics that controls the shape of 
these corrections is by definition unknown, it seems impossible to establish a generic result with regards to the shape 
of the modified power spectrum. Here we will consider a conservative starting point, namely that the perturbations 
begin in the expanding phase in their instantaneous minimum energy state (as we have seen before, this is not the 
most general case). 

We have seen before that, if trans-Planckian physics is adiabatic, then the inflationary initial conditions are not 
modified. This means that the sub-Hubblc Mukhanov-Sasaki variable is still given by 

V2k m pl V2T m Pi 

with a;, = 1 and fit. = 0. If, however, the physical conditions that prevailed in the trans-Planckian regime were 
non-adiabatic, then there was particle production and, as a result, (3k ^ 0. This immediately implies that the power 
spectrum (which is proportional to the square modulus of Vk) contains super- imposed oscillations. Therefore, even 
without knowing in detail the trans-Planckian physics, it is possible to establish a generic prediction for the shape of 
the modified power spectrum. This generic prediction was made for the first time in Rcfs. [19]. 

However, clearly, if we want to make more precise predictions and, for instance, compute the amplitude, the 
frequency or the phase of those oscillations, we need to make further assumptions. The most general approach 
consists in parameterizing the deviations from the standard situation in the initial conditions. Let M c be the energy 
scale above which some new physics is operating. This scale could be the Planck scale or the string scale for instance. 
A Fourier mode emerges from the regime where the new effects are relevant when its physical wavelength equals the 
new length scale introduced in the problem, namely 

m = T a(v)=ic = Wc , (ID 

where k is the comoving wavenumbcr. It is important to notice that one can have A <C £ c and, at the same time, 
H <C m pl where H is the Hubble parameter during inflation. In other words, the Fourier mode wavelength can be 
much smaller than the new fundamental scale in a regime where spacetime can still be described by a classical FLRW 
background. The initial conformal time satisfying Eq. (11) depends on the scale k and, for this reason, it will be 
denoted in the following by r\k- At this time, if the trans-Planckian regime was non-adiabatic, then we expect the 
initial conditions to be modified compared to the standard case. Technically, this is expressed by 

, , a k + /3 k 40F 
Vk(Vk) = , (12) 

VM^mJ to pi 

i l \ j / ^(7]fc)40F(q fc -/3fc) 
v k {Vk) = \ — ^ , (13) 



where ak and /3 k are two complex numbers (already introduced before) satisfying \ak\ 2 — |/3fc| 2 = 1- These numbers 
completely characterize the influence of the new physics on the initial conditions. Again, a priori, a complete calcu- 
lation of ak and /3k requires the knowledge of the physics beyond the scale M c . However, if the ratio H/M c goes to 
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zero, then we expect to recover the standard Bunch-Davies situation for which a& = 1 and = 0. Therefore, if the 
expressions of a*, and j3k & re perturbative in H/M c , then, without loss of generality, one can write 
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where x and y are two numbers characterizing the perturbative expansion. The modified power spectrum can be then 
determined perturbatively in terms of the free parameters x and y. Notice that, strictly speaking, x and y can be 
scale-dependent. Here, for simplicity, and in order not to introduce arbitrary functions in the problem, we assume 
that they are roughly scale- independent over the range of wavenumbers relevant to the CMB. If this is not the case, 
then one looses the ability to parameterize the modified power spectrum in a simple way and one would really need a 
complete theory of the trans-Planckian effects to establish the shape of the corrections. It is also important to recall 
that afc and (3k must satisfy the Wronskian condition. At leading order in H/M c , this simply amounts to y + y* = 0. 

We are now in a position where the modified power spectrum can be derived. Following Ref. [32] , it can be written 
as (of course, there is a similar formula for tensor perturbations that we do not show here) 
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In this expression, t\ and e 2 are the two first slow-roll, or horizon-flow, parameters. The scale fc* is the pivot scale 
and ao is the scale factor evaluated at a time r/o during inflation (the Hubble parameter which appears in the overall 
normalization is also evaluated at the same time). This time is arbitrary and does not depend on the scale fc. For 
convenience, we choose fc*/ao = M c . The quantity C is given by C = j E + In 2 — 2 where 7 E is the Euler constant. 
Finally, ip is the phase of the complex number x, that is to say x = \x\c lv . It is interesting to notice that, at this 
order, the number y does not appear in the expression for k?P^. 

Let us now comment on the power spectrum itself. We see it has the expected shape. There is the usual slow-roll 
(almost scale-invariant) part, with the overall normalization given by the square of the Hubble parameter (during 
inflation), measured in Planck units, and divided by the first slow- roll parameter. This slow-roll part is corrected by 
two oscillatory terms which represent the super-imposed oscillations. It is interesting to notice that these oscillations 
are logarithmic in fc. This is due to the fact that the initial conditions have been chosen on the surface M c = 
constant [68, 69]. If, on the contrary, the initial conditions were chosen on a surface of constant time, then the 
oscillations would be linear. This is for instance what happens in bouncing models. The amplitude of the oscillations 
is, roughly speaking (taking |x| ~ 1), given by \x\H/M c while the frequency is proportional to (7J/M C ) _1 . Therefore, 
if the initial state is the instantaneous Minkowski vacuum, the amplitude is just inversely proportional to the frequency. 
This means that high frequency modes necessarily have a small amplitude. But, of course, there is a priori no reason 
to assume that the frequency and the amplitude are fully correlated. In addition, as we are going to see, postulating 
\x\ = 1 would prevent us to find the best fit. 

On general grounds, we expect the ratio H/M c to be a small number. Indeed, we know from the COBE normal- 



ization that H < 10 5 r 



We do not know the scale M„ but reasonable candidates are the Planck scale itself or 



the string scale M s ~ 1CP 1 — 10~ 3 to p1 . This means that H/M c is at most of order ~ 0.01. Therefore, if |x| = 1, the 
amplitude of the trans-Planckian corrections is very small and it will be difficult to detect an effect in the CMB, even 
with an experiment like the Planck satellite. 

On the other hand, if |x| ^ 1, one can easily reach a level which would be detectable. In this case, we nevertheless 
meet another issue, namely the backreaction problem [70, 71]. If \x\ ^ 1, we start from a non- vacuum state, a natural 
assumption if physics beyond the scale M c is non-adiabatic. But a non- vacuum state means that particles are initially 
present. Since these particles carry energy density, there is now the danger that this energy density be comparable to 
or larger than that of the background, (which is ~ H 2 m? >] ), in which case the energy in the particles would prevent 
the onset of inflation. Clearly, the larger \x\ is, the larger the trans-Planckian corrections become, and the more severe 
is the danger to have a backreaction problem. In fact, it is easy to estimate the regime where there is no backreaction 
problem. The energy density carried by the "trans-Planckian particles" is 
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FIG. 6: Trans-Planckian power spectra given by Eq. (16). The blue line corresponds to a vanilla model with |a;| = and 
ei = 1/(2AAT,), €2 = 1/ANt with AiV* ~ 50 as predicted for the m 2 (f> 2 inflationary model. The red line corresponds to a 
model with the same values for the slow-roll parameters and H/M c ~ 0.002, \x\ ~ 50, ip — 3. Finally, the green line represents 
a model with H/M c ~ 0.001, ip = 2 and the same values for the other parameters. 



Therefore, the condition p uv < i? 2 m p[ reduces to \x\ < m pl /M c or, using the COBE normalization H ~ m pl 10 4 -y/eT, 



As a consequence, one can have \x\H/M c of order one and easily satisfy the above inequality (in particular, in 
inflationary models where the gravity waves contribution is extremely small). Of course, this does not mean that one 
should see trans-Planckian effects in the CMB. It could very well turn out that |x| is not large enough to compensate 
the smallness of the ratio H/M c . In that case, the trans-Planckian effects arc there but it seems difficult to imagine 
how the small oscillations could be measured one day. 

Another remark is in order here. In the above calculation, we have used the value of the Hubble parameter during 
inflation. As noticed in Ref. [72], if instead one had used its present day value, the no-backreaction condition would 
have become almost impossible to satisfy. Using the Hubble parameter during inflation assumes that the trans- 
Planckian production of particles stops after the end of inflation. It is easy to find a situation where this happens. 
Let us for instance consider the case where the trans-Planckian corrections arise from a modified dispersion relation 
(see the previous sections). A typical situation is represented in Fig. 7 (the model and the argument presented here 
were also studied in Ref. [73]). The dispersion relation is almost linear for k v \ ] y a /m c <C 1 as it should be for obvious 
phcnomcnological reasons, and deviates from linearity for kp\ lys /m c I. In the regimes where w p hys ^ H, the 
adiabatic approximation is violated and particle production takes place. This is of course the case in the regime 
^phys during inflation, this situation corresponding to the case where the wavelength of a Fourier mode is larger 
than the Hubble radius. But, as can be noticed in Fig. 7, it is also the case when fc p h ys G [k-, fc+] during inflation. 
In this regime, adiabaticity is violated and the resulting inflationary power spectrum is modified (and, typically, as 
argued before, acquires super- imposed oscillations). But the Hubble parameter evolves (decreases) during inflation 
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FIG. 7: Example of a non-linear dispersion relation leading to a modified power spectrum. In the regime fc p h ys <C M c , the 
dispersion relation is approximatively linear, as appropriate for standard massless excitations, while in the regime fc p h ys 3> M c , 
it strongly deviates from linearity due to trans-Planckian effects. In the regions where u> p h ys < H, the adiabatic approximation 
is violated and trans-Planckian particles production occurs. The horizontal dashed red lines represent the value of the Hubble 
parameter at different times during the cosmic evolution. In this model, particles production happens during inflation for 
modes such that fc p h ys < fco and fc p h ys £ [k-, k+] and, after inflation, only for modes fc p h ys < fco, this latter regime taking place 
in the linear part of the dispersion relation and corresponding to the standard amplification of cosmological perturbations on 
super-Hubble scales. 



and after. As a consequence, the horizontal red dashed line in Fig. 7 representing H/M c goes down in this diagram. 
At some point, it passes below the minimum of the dispersion relation and particle production stops. In particular, 
this is the case for the line representing the current Hubble parameter. In such a situation, it is clear that -ffinf (and 
not Hq) should be used in order to evaluate p vv . 

Moreover, we want p vv <C -H^mj^ because the energy density of the "trans-Planckian particles" could spoil inflation. 
But, of course, this rests on the prejudice that the corresponding equation of state is different fromp/p = —1, typically 
that of radiation p/p = 1/3. When the dispersion relation is modified, as shown in Refs. [72], the equation of state 
also acquires corrections. And, as was demonstrated in Rcf. [72], in the case of the dispersion relation in Fig. 7, the 
modified equation of state is precisely very close to that of the vacuum. Therefore, even if one has a backreaction 
problem, it is not obvious that this will prevent inflation to proceed. 

Before turning to the observational constraints, it is also interesting to study how the oscillations in the power 
spectrum are transferred to the multipole moments. In principle, this calculation can only been done numerically. 
There is, however, a regime where one can obtain an approximate analytical result. For small £ (that is to say large 
angular scales), in the limit t\M c jH 3> 1 and for ei — — 2ei (or n s = 1), the multipole moments can be expressed 
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FIG. 8: Multipole moments in presence of super-imposed trans-Planckian oscillations, courtesy of Ref. [86]. As discussed in 
the text, the high-frequency super-imposed oscillations are damped on large scales and appear at the rise of the first Doppler 
peak. 
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where r\ ss is the distance to the last scattering surface. The first conclusion that can be drawn from the above equation 
is that super-imposed oscillations in fc-space are indeed transferred to the €-space, as the presence of the trigonometric 
function shows. This type of feature is not wiped out by the transfer function because the oscillations are present 
everywhere in fc-space. Moreover, we see that logarithmic oscillations in fc-space give rise to linear oscillations in real 
space since the argument of the cosine function scales as ~ tt£. In fact, this last property turns out the be an artefact 
of the approximations used here. Since we are in a situation where the frequency of the oscillations in momentum 
space goes to infinity, one finds oscillations in £-space with frequency 7r which is just the maximum frequency possible 
given that the values of £ are discrete. 

Finally, the amplitude of these oscillations is damped by a factor (eiM c /H) 5 / 2 for very small £, this effect being 
compensated by £(£+ 1) for larger £. However, since the above formula breaks down in the latter regime, it is difficult 
to reach definite conclusions. One can nevertheless say that we expect the super-imposed oscillations to be absent 
at very large scales and to appear only at relatively small scales. These considerations are confirmed in Fig. 8 where 
multipole moments obtained from a trans-Planckian power spectrum have been displayed. In particular, we see that 
the super-imposed oscillations only appear at the rise of the first peak but are damped at very large scales. 

Having studied how the trans-Planckian effects affect the CMB anisotropy, let us now turn to the question of their 
possible detection in astrophysical data. 




FIG. 9: Multipole moments measured by the WMAP experiment. The calculation of the Ci depend on the initial power 
spectrum and, therefore, the corresponding data can be used to constrain trans-Planckian physics. 



IV. OBSERVATIONAL CONSTRAINTS 



A. The Power Spectrum 



We are now ready to discuss the observational constraints on trans-Planckian physics. It is conventional to assert 
that quantum gravity can be probed by studying the propagation of high energy cosmic rays. However, this is not 
the only astrophysical method which allows us to constrain trans-Planckian effects. This can also been done with the 
help of the CMB anisotropies that have been precisely measured by the WMAP satellite [75, 76] (see Fig. 9). Here, 
we briefly review the main techniques that have been used to address this question [74, 77, 78] (for related works in 
the context of inflation without trans-Planckian effects, see Refs. [79-83]). The exploration of the parameter space 
corresponding to Eq. (16) is usually done using Monte Carlo techniques. Each model is computed using a modified 
version of the CAMB code [84] and the likelihood is estimated using the COSMOMC code [85]. 

A first problem met in this kind of analysis is that, in presence of super-imposed oscillations, the calculation of 
one set of multipole moments becomes very time consuming (a few minutes instead of a few seconds). As a 
consequence, it is not possible to explore the full parameter space and the analysis must be restricted to the "fast 
parameter space" . This means that the bare cosmological parameters are fixed to their best fit values and that the 
Monte Carlo exploration is only performed on the "primordial parameters", namely the overall normalization of the 
spectra, the slow-roll parameters and the parameters describing the oscillations. 

The parameters describing the super-imposed oscillations are the amplitude, the frequency and the phase. The 
parameter describing the amplitude is, as discussed above, \x\H/M c and it is convenient to take a uniform prior on 
this parameter in the range [0, 0.45]. In order to sample the frequency, one can choose a uniform prior in [1, 2.6] for 
the parameter log(eiM c /H). There is no need to emphasize again the important role played by the fact that \x\ =/= 1. 
Taking |a;| = 1 is not only poorly theoretically justified but, from the data analysis point of view, would render the 
amplitude and the frequency completely degenerate, a property that can affect a lot the result of the Monte Carlo 
exploration. Finally, the phase is described by the parameter tp = 2M C (1 + (-\)/H + ip and one assumes a uniform 
prior in [0, 2tt]. 

The converged posteriors, marginalised posteriors (solid black line) and mean likelihood (purple shaded bars), for 
the amplitude, frequency and the phase, have been plotted in Fig. 10 given the seven years WMAP data [86]. The 
marginalised posterior is a quantity which is sensitive not only to the absolute value of the likelihood function but 
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FIG. 10: Posterior distributions for the three trans-Planckian parameters describing the frequency, amplitude and phase of 
the super-imposed oscillations, courtesy of Ref. [86]. The parameter <tq stands for H/M c . The solid black lines represent the 
marginalized distributions while the purple shaded bars represent the mean likelihood. 



also to the volume occupied in parameter space. Therefore, if a set of parameters leads to a very good fit but is 
such that it requires a precise fine-tuning of those parameters (that is say, if one slightly detunes the parameters, 
then the likelihood dramatically drops), then its marginalized probability can be small. On the other hand, the mean 
likelihood, as its name indicates, is only sensitive to the absolute value of the likelihood function and is independent 
of volume effects in the parameter space which is being considered. 

Let us now describe the results displayed in Fig. 10. The amplitude of the marginalised probability is peaked over 
a vanishing value which indicates that the vanilla slow-roll power spectrum (with no trans-Planckian oscillations) is 
still the favored model from a Bayesian point of view. One finds that \x\H/M c < 0.26 at 2er confidence level. It is, 
however, also interesting to remark that, as opposed to the marginalised probability, the mean likelihood is peaked 
over a non-vanishing value. This means that the best fit is actually a model with trans-Planckian oscillations. One 
finds that A% 2 ~ — 11 compared to the vanilla slow- roll model. The discrepency between the marginalized probability 
and the mean likelihood is interpreted, in agreement with the above discussion, as an indication that volume effects 
play an important role. In other words, if one goes away from the best fit model, the likelihood function quickly 
decreases which means that the best fit only occupies a small volume in parameter space. This analysis is confirmed 
by the frequency posteriors which exhibit a series of peaks at particular frequencies. The fact that those peaks are 
very narrow indicates that it is necessary to fine tune the frequency parameter in order to find the best fit, which is 
totally consistent with the fact that the best fit model occupies a small volume and, hence, that the marginalised and 
mean likelihood posteriors differ. Let us add that, in Ref. [87], it was shown that, if the data actually contain super- 
imposed oscillations, then the likelihood function should precisely exhibit an oscillatory structure similar to the one 
which seems to emerge from the present analysis (this structure is even more clearly visible on the 3D representation 
of the posteriors plotted in Fig. 11). This is an intriguing fact. Finally, the frequency remains basically unconstrained 
although ip ~ 6 seems slightly favoured but, in any case, not in a statiscally significant way. 

We have emphasized previously that the trans-Planckian effects lead to super-imposed oscillations that are loga- 
rithmic in k. It is therefore interesting to test whether this particular scale dependence is favored by the data. This 
was done in Ref. [79] for the WMAP3 data. The idea is to postulate the following power spectrum 
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such that, in the limit p — > 0, one recovers logarithmic oscillations. Assuming a uniform prior on logp in the range 
[—5, 0.48], one obtains at 2a confidence level that p < 0.68. In this sense, one can indeed argue that the logarithmic 
structure is special. 

The above discussion should however be toned down given the fact that the best fit suffers from a severe backreaction 
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FIG. 11: Three-dimensional marginalized probability distributions in the amplitude- frequency space (left panel) and in the 
amplitude-phase space (right panel), courtesy of Ref. [86]. 

problem. Indeed, the best fit is obtained for the following values: \x\H/M c ~ 0.13, \og(eiM c /H) ~ 2 and log(ei) ~ 
—3.1 (of course, this last result does not imply a detection of primordial gravitational waves since we only explore 
the fast parameter space, see the discussion before). These numbers imply that H/M c ~ 10~ 5 and \x\ ~ 10 4 . It 
is interesting to notice that the ratio H/M c has a very natural value and that it was indeed very useful to perform 
the analysis without the assumption \x\ = 1. However, Eq. (18) now reads \x\H/M c < 1CT 4 - 5 and the best fit 
clearly violates this inequality by three or four orders of magnitude. Of course, this does not invalidate the statistical 
analysis presented here, it just questions the interpretation of the best fit in terms of super-imposed trans-Planckian 
oscillations. 

So, in conclusion, what is the observational status of the trans-Planckian oscillations? Clearly, the best model 
remains the slow-roll vanilla model and there is no detection, at a statistically significant level, of super-imposed 
oscilations. The Bayesian evidence for this model has never been computed (it was computed for the first time for 
different slow-roll models only recently in Ref. [83] ) but it seems pretty clear that its estimation would only reinforce 
this conclusion. Moreover, we do not see any clear trend towards a detection as more and more accurate data set are 
released. Indeed, for the first year WMAP data, we had \x\H/M c < 0.11, for the three year WMAP data we obtained 
\x\H/M c < 0.38, which is a "better" result, but the result was \x\H/AI c < 0.26 for the seven year data. There are, 
however, a list of intriguing hints (the best fit is given by a model with super-imposed oscillations, presence of narrow 
peaks in the likelihood function) that seems to indicate that, maybe, non trivial features are actually present in the 
data. Of course, if this is the case, it remains to be seen whether there are of primordial origin and whether they 
can be explained by trans-Planckian physics. But we are of the opinion that these hints are a sufficient motivation to 
carry out the analysis again when the Planck data become available. 

B. Non Gaussianities 

Let us now study how trans-Planckian effects can modify the conventional predictions for non-Gaussianities. As is 
now well-established, the level of non-Gaussianity is small in single field inflationary models (with a canonical kinetic 
term) if the initial state is the Bunch-Davies state. Roughly speaking, as shown for the first time in Refs. [88-92] 
and then further elaborated in Refs. [93], the / NL parameter (sec the definition below) is of the order of the slow-roll 
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parameters 5 . This level is too small to be detected even with the Planck satellite (the expected level of detection 
i s I/nlI ^ 5). We have seen, however, that trans-Planckian physics, if non-adiabatic, can place the cosmological 
perturbations in a non- vacuum state for which the above-mentioned result no longer holds. It is therefore interesting 
to estimate the level of non Gaussianities if trans-Planckian effects are present in the early Universe. 

The first calculation in this direction was performed in Refs. [94-96]. It was done for a non- vacuum state with 
a "preferred direction" in momentum space such that the three-point correlation function does not vanish even in 
the free theory. In this case, the level of non-Gaussianity was found to be extremely small and beyond the level 
of detectability. We are, however, more interested in computing the non-Gaussianity for a state characterized by 
Eqs. (14) and (15) which generically describe the trans-Planckian modifications. In this case, the free theory is still 
Gaussian and, therefore, the corresponding three-point correlation vanishes. It is only when interactions are taken 
into account (that is to say beyond the Gaussian approximation) that non-Gaussianity can be produced. This problem 
has recently been considered by various authors, see Refs. [97-101], and, in the following, we turn to this question. 

A Gaussian distribution has a vanishing three-point correlation function and any detection of a non-vanishing 
signal would therefore rule out Gaussian perturbations. For this reason, it is interesting to calculate the three point 
correlation function of the curvature perturbation defined by 

<C(*»K(**Kfo,*)> = J ^ J^- 2 / ^<Oa(„)0 M Wa.W>e«t w +»^»-. (21) 

In the above expression represents the Fourier transform of £ (and kl = k\ etc.). Let us recall that the curvature 
perturbation ( is related to the Mukhanov-Sasaki variable v introduced before by C = v / (y/2M P} a,y/ei) . In the 
framework of the theory of cosmological perturbations of quantum-mechanical origin, it is an operator and it can be 
written as 
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Since one has Cfe = (ahfk + a kfk), the power spectrum can be expressed as [see also Eq. (2)] 
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where the quantity k 3 P<^ was also already introduced before, see for instance Eq. (16). The most general form of Vk 
is given by 
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Several comments are in order here. Firstly, the coefficients ak and /3k are of course the same coefficients used to 
evaluate the trans-Planckian power spectrum. The Bunch-Davies vacuum corresponds to the choice ak = 1 and 
Pk = 0. Here, we are interested in a situation where /3fc 7^ and a generic parameterization of the deviations from 
the Bunch-Davies state due to trans-Planckian effects has been presented in Eqs. (14) and (15). Secondly, one can 
check that the above expression behaves according to Eq. (10) in the sub-Hubble regime (in fact, we have slightly 
changed the conventions - no factor Ay/iv/m pl appears anymore - in order to work with notations that are standard 
in the calculations of non-Gaussianities). Thirdly, Eq. (24) corresponds to a de Sitter background. This means that 
we neglect the corrections due to the slow-roll parameters in the mode function (but we took them into account when 
we computed the power spectrum). These corrections would lead to sub-leading effects in the calculation of the three 
point correlation function (for a calculation, see for instance Refs. [102, 103]). Finally, from Eq. (24), it is easy to 
establish the expression of fk that will be used later on. One finds 

t k = iakI ^ (l + ik v )c- 1 ^ ^$5=(1 - i^tfe**". (25) 



5 Although this is very often not properly acknowledged in the literature, the fact that the non-Gaussianity is small for the class of models 
mentioned above was established in Refs. [88-91] much before the publication of Ref. [93]. The main contribution of Ref. [93] has been 
to calculate the momentum dependence of the bi-spectrum exactly and to introduce efficient techniques which have then permitted the 
calculation of non-Gaussianities for other, more complicated, classes of inflationary scenarios. 
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We will also need the derivative of the mode function. One obtains 
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where we have neglected terms suppressed by the slow-roll parameters. 

As we discussed before, the three point correlation function can only be non- vanishing if interactions are taken into 
account. In this case, the three point correlation function in momentum space (or bi-spectrum) can be expressed 
as [93, 104] 



{Cki(v) 6=2(77) Cfc3(?7)) = -if dra(r) ([Chi(v) Ck2(v) Ch3(v), H iat (r)} 



(27) 



with fli n t being the interaction Hamiltonian while r)i n i is the time at which the initial conditions are imposed on the 
modes when they are well inside the Hubble radius. On the other hand, rj c denotes the time when inflation ends. In 
practice, we always take rj e — > 0. The choice of ?7ini is more tricky in the trans-Planckian context and will be discussed 
in more detail in the following. In order to calculate the interaction Hamiltonian, one must compute the action for 
cosmological perturbations at cubic order (the action of the free Gaussian theory is quadratic). A now standard 
calculation gives [93, 102-104] (a dot means derivative with respect to cosmic time) 
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where 8C2/5C denotes the variation of the second order action with respect to C tl and is given by 
and the quantities A and x are defined by 

a 2 j> 2 



A 



- l C = a 2 e 1 C, x = d~ z A. 



2NP H 2 

PI 

The term J : (SC2/SC) introduced before refers to the following complicated expression 



(28) 
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(31) 

Moreover, the terms which involves SC2/SC can be removed by a suitable field redefinition of C of the following 
form [93, 102-104]: 



C -> Cn + F(Cn). 



(32) 



With such a redefinition, the interaction Hamiltonian corresponding to the action 53(C), an d to be used in Eq. (27), 
can be written in terms of the conformal time coordinate 77 (a prime means a derivative with respect to conformal 
time) as 



flint (V) = d " X 



4CC' 2 +aelC(dC) 2 -2e 1 C (5*0 



+ \ 61 e ' 2 eC + Y~a {dX) {d * X} ^ + Ta ( ^ C) {9xf 



(33) 
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The three first terms are second order in the slow-roll parameters while the three last ones are third order. Despite 
the appearence of the slow-roll parameters, this expression is exact at third order in the perturbations. 

Having determined the interaction Hamiltonian, one is now in a position where non-Gaussianities can be calculated. 
For convenience, we redefine the bi-spectrum according to 



(Cfcl0?e) Ck2{Ve) Cfc3(%)> = 



(27T) 



9/2 



G(kl, k2, fc3) <5 (3) (fel + k2 + kS) . 



(34) 



where the delta function ensures momentum conservation. Then, the quantity Q(kl, fe2, fc3) can be written as 



£?(fcl,fc2,fc3) = 



c=\ 



fki i.Ve) fk 2 (Ve) fk 3 (Vc)G c (fel, fe2, fe3) 



-fti foe) f k2 (Jh) fl (Ve)G* c (fel, fe2, fe3) 



(35) 



The quantities Q c (fel, k2, k3) with C = (1,6) correspond to the six terms in the interaction Hamiltonian (33), and 
are given by [93] 

(?i(fel, fe2, fe3) = 2i I dr a 2 e\ (/^ f£ 2 fuz + two permutations) , 
J mm 

£? 2 (fel,fe2,fe3) - 
g 3 (fel,fe2,fe3) 



2i / dra e 1 f^ 2 fk3 (fel ' fe2 + two permutations) 

rim 



2i I dr a 2 e\ 



/* rf* rl* 
fel /fc2 Jk3 



fel • fe2 

k 2 



£? 4 (fcl,fe2,fc3) 
£ 5 (fcl,fc2,fc3) 

£ 6 (fel,fc2,fc3) 



+five permutations 

f»7e 

,2 , J 



(36) 



dr a ei e 2 (/^ /fc 3 + two permutations) . 



dr a 2 ef 



/* rl* /V* 
fel Jfc2 /fc3 



/fel ■ fc2 



+five permutations 



dr a 2 e 3 



/* W* j 1 /* 
fel ifc2 Jfc3 



(fe2 • fe3) 



■ two permutations 



Actually, an additional, seventh term arises due to the field redefinition (32), and its contribution to G{kl, fc2, fe3) is 
found to be 



g 7 (kl, fc2, fe3) = y (|/ fc2 | 2 |/fc 3 | 2 + two permutations) 



(37) 



The other terms in Eq. (31) do not contribute because they all contain a derivative (time derivative and/or space 
derivative) and, at the end of inflation, on large super Hubble scales, C is constant. 

The last thing which remains to be done is to define the / NL parameter which is conventionaly used in the literature 
to quantify the level of non-Gaussianity. It reads 



(CfeiCfe2Cfe3> = (2tt) 4 + 3 / 2 (-A /nl7 ^E 



1.3 



S(kl + k2 + k3), 



(38) 



where we have neglected the slight scale dependence of V s . This definition is introduced to match the formula 

C(ry, x) — £ G — 3/ NL (£ G ) /5 which is used to extract observational bounds on the parameter / NL (the coefficient 3/5 
comes from the fact that, during a matter dominated era, the curvature perturbation and the Bardcen potential are 
related by C = 5$/3). 
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^From the above expression, the full bi-spectrum can be calculated. For a slow-roll model with canonical kinetic 
term, the dominant terms are Q\, Q2, S3 and Q7. This may differ for other type of models. For instance, if slow-roll 
is violated during a few e-folds, then the main contribution comes from the term O4, see Ref. [105-108]. If the kinetic 
term is not standard, then a new vertex cx £ 3 appears, see Ref. [102-104]. Here, we do not give the full bi-spectrum for 
a slow-roll model since this is not the main goal of this review article. It can be found in Refs. [93, 102-104] (Notice 
that it was recently shown that this bi-spectrum is not modified by re/pre- heating effects in single field inflation, 
see Ref. [109]. Therefore, it actually corresponds to what we should see in the CMB). We just quote the so-called 
consistency relation for the squeezed configuration k\ ~ fc 2 3> k 3 which reads 

/S = ^K-1), (39) 

and expresses nicely the fact that the non-Gaussianity is indeed of the order of the slow-roll parameter (it is also 
important to recall that the validity of this formula is in fact broader, see for instance Ref. [110]; in the case of ultra 
solw-roll inflation, it is violated but the corresponding solution is in fact instable [111]). 

Having briefly reviewed the standard situation, let us now turn to the evaluation of the bispectrum in a situation 
where (3k 7^ 0. In principle, one could just restart from Eq. (28). However, it was noticed in Ref. [93] that the first 
three terms in this equation (that is the say the three dominant terms, second order in the slow- roll parameters) can 
be reduced to 

S3 [C] = Ml J dt d 3 x Ac?HelQ' 2 d- 2 £ , (40) 

up to higher order terms in the slow-roll parameters and terms involving 8C2/ &(,■ This can be achieved by performing 
various integrations by part. The calculation is then reduced to the calculation of a single vertex. This means that 
the expression (35) is still valid but that the three vertices Q\, Q2 and Q3 are now replaced by a single vertex G123 
given by 



'123 




' k dra 3 Hel f&f&f&. (41) 



In addition, since we have seen that G4, G5 and Qq arc in fact sub-dominant in the slow- roll approximation, we only 
have to calculate the above vertex. Using the expressions (25) and (26) for the mode function and its derivative, one 
arrives at 

(cm cm cm) = (2-r 3/2 <s (fci + k2 + fc 3 ) y^-jpM. i] 

° m P1 e l [[ l=1 Kj \ i=1 K i I 



aki - Phi) («fc2 " /3fc2) (ak3 ^ Pm) / dr (a* kl c lklT - fi* kl c 

J mm 

^ 2 e- lfe2T ) K 3 e* fe3r - PLe- lk3T ) c.c 



(42) 



We can check that, if /3k = in the above equation, together with the vertex £7, one recovers the consistency 
relation (39). If /3fc 7^ however, interesting new effects appears. For instance, typically, the bi-spectrum contains the 
following term 

(CM CM CM) ^ «li PI2 a* k 3 r dre^- k ^ 



(43) 

mm 



i(k 1 -k 2 +k3) 



i{ki - k 2 + k 3 ) 



We see that, in the squeezed limit k\ = k2 ^> ^3, there is an enhancement of the signal by a factor ki/k$ ^> 1 [100]. 
One can therefore hope to constrain excited states (and hence trans-Planckian physics) with the help of the future 
non-Gaussianity measurements. 
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In fact, it is straightforward to calculate the full bi-spectrum for the generic parameterization given by Eqs. (14) 
and (15). One obtains 

(Cfcifa,) CM CM) = (2n)-^6 (fcl + fc2 + fc3) ^ fei) 

H H 

- 3|s| — cos V3 + 3|a;| — cos (k T i] ini + (p) 
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H k 1 

■cosv? + 3|a;|— — ^ cos [(-fci + fc 2 + fc 3 ) Jftoi -<p]>, (44) 

M c -fci + fc 2 + fc 3 J 

where we have defined fc T = fci + fc 2 + k 3 . This expression is valid at leading order in When |x| — > 0, one recovers 
the usual result. Let us notice that, in this case ?yi n i should be rotated to the complex plane and sent to infinity. With 
this usual procedure which singles out the vacuum state of the theory, the term cos (fc T 77; n ;) also vanishes. The above 
expression is quite complicated but we see that the bi-spectrum is enhanced when fcj = X)i=i — ^fc; vanishes [98]. 
In order to have a clearer view of the result, it is interesting to consider the squeezed limit again. In the limit 
fci = fc 2 = fc 3> /C3, the above equation takes the form 

(Cki(Vc)Ck2(ri c )Ck 3 (r]c)) sq = (27r)- 3 / 2 <5(fcl + k2 + fc3) 



H k 

ei + 6ei|x|— -— cos (k 3 r) irli - tp) 
M c fc 3 



(45) 



This simple expression captures the main features of the result. We see that the corrections are proportional to 
\x\H/M c . We also see the enhancement in the squeezed limit, proportional to fc/fc3 3> 1. This brings us to the 
question of the initial time ?yi n i. We have seen that, in the minimal approach, the modes are created at different times, 
when their wavelength equals the new fundamental scale M" 1 , see Eq. (11). In an inflationary background, this implies 
that ?7i n i = r]k — M /(kH), that is to say the initial time becomes scale dependent. However, when one computes the 
various vertices that contribute to the bi-spectrum, one has to integrate the product of three mode functions with 
different wavenumbers fcl, fc2 and fc3 from rj[ n [ to zero, see the previous expression for Q123 for instance. Then comes 
the question of which wavenumber should be considered in the expression of ryi n ; = rjk — M c /{kH)7 Let us notice, 
however, that if the scaling ?7i n i cx M c /(kH) is correct (whatever the precise definition of fc is), then the frequency 
of the oscillatory term in the bi-spectrum becomes cx M c /H in a way which is very similar to the super- imposed 
oscillations found in the power spectrum, see Eq. (16). 

Finally, let us remark that once the bi-spectrum has been determined in momentum space, one still needs to perform 
the 2D projection in order to obtain the / NL seen in the sky. This is a highly non trivial procedure [98, 100, 112]. 
According to these references this could lead to an observable / NL which, therefore, could constrain trans-Planckian 
physics. 



V. BROADER CHALLENGES 



A. UV and IR Cutoff Issues 



Let us now return to some more general issues related to trans-Planckian physics in cosmology. The trans-Planckian 
problem for cosmological perturbations is related to basic challenges to the applicability of effective field theory in an 
expanding background 6 . 

In free quantum field theory, the starting point is canonical quantization of the fields. In an expanding background 
(see [33] for an overview of quantum field theory in curved space-times), the basic modes which are quantized are 



Similar issues arise in applications to black hole backgrounds and to non-gravitational external field background. 
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the plane wave modes which have constant wavelength in comoving coordinates and whose physical wavelength is 
hence increasing. Even in Minkowski space-time, quantum field theory requires regularization and renormalization 
in order to obtain finite answers for observables. In the context of gravity, pure quantum field theory must break 
down on wavelengths smaller than the Planck scale since waves with such large wavelengths would trigger the collapse 
of the background space-time into a gas of black holes. This problem is intimately related to the cosmological 
constant problem which arises when considering quantum field theory in gravity. Hence, in the presence of gravity, 
the ultraviolet cutoff is not just a computational tool, but it as an issue with definite physical significance. 

The problem in an expanding background is that the ultraviolet cutoff is expected to correspond to a fixed scale 
in physical coordinates, whereas the basic modes of a quantum field have constant comoving wavelength. Thus, as 
already discussed in [64], a fixed physical UV cutoff implies that the Hilbert space of the quantum field must be 
time-dependent. As the universe expands, more and more new modes are required to describe the physics. 

The above UV problem arises in any cosmological space-time. What is special to inflationary cosmology (as opposed 
to the alternatives which we have mentioned in Section 2) is that these new modes are inflated to a wavelength which 
at the current time is observable. 

For masslcss fields, there are also infrared (IR) divergences. Since gravitational waves are masslcss, potential 
infrared divergences must be addressed in cosmology. 

In cosmology, there is a natural scale which separates UV and IR modes - the Hubble scale £f _1 (i), where H(t) is 
the expansion rate of space at time t. On sub-Hubble scales matter fluctuations oscillate as they do in flat space-time 
7 . However, on super-Hubble scales the oscillations are frozen out and plane wave fluctuation modes become standing 
waves whose amplitude evolves in time as determined by the gravitational background. Furthermore, a local observer 
(local in space and time) has a range of view limited by the Hubble radius. Hence, the Hubble radius determines the 
separation between UV and IR modes. 

The increasing physical wavelength of the basic modes of a quantum field also leads to an infrared problem: In the 
case of an accelerating universe such as in inflationary cosmology, the phase space of infrared modes increases. Even 
in the alternatives to inflation which were mentioned in Section 2 this problem arises. In fact, the increase of the 
phase space of infrared modes is a necessary condition to have a causal theory of structure formation since we need 
to make sure that scales which are being observed in cosmology today (and which were in the infrared sea, i.e. had 
a wavelength greater than the Hubble radius for most of the late time universe) start out with a wavelength which 
is sub-Hubblc. In inflationary cosmology, the phase space of infrared modes increases during the inflationary phase, 
in string gas cosmology it grows at the end of the Hagcdorn phase, and in a bouncing cosmology it grows during the 
contracting phase. 

Infrared divergences can be "cured" by imposing an infrared cutoff. Such an infared cutoff can be well justfied on 
physical grounds: it corresponds to setting to zero the contribution of modes which are super-Hubble at the initial 
time, and whose values would depend on dynamics before that time. In an expanding universe, this cutoff must be 
fixed in comoving coordinates, in contrast to the UV cutoff which is fixed in physical coordinates (for a discussion of 
this point see e.g. [113]). 

Since the phase space of IR modes is increasing relative to the phase space of UV modes, the magnitude of the 
infrared effects will be time-dependent. A possible consequence of this will be discussed in the following subsection. 

As mentioned above, the growth of the sea of IR modes occurs in any causal theory of structure formation. What is 
special in inflationary cosmology is that the sea becomes populated by modes which initially were on trans-Planckian 
scales. Thus, in inflationary cosmology the trans-Planckian problem for cosmological fluctuations and the infrared 
issues become coupled, which they are not in string gas cosmology or in a matter bounce scenario. 

The effects of infrared modes on the spectrum of cosmological perturbations in inflationary cosmology has been 
studied recently in a large number of works [114]. A first point to consider is that the effects of IR modes is different 
in an inflationary universe (a cosmology with a finite duration dc Sitter phase) compared to what happens in exact 
de Sitter space. As recently studied in [115], the functional form of the IR terms is different in the two cases, the 
extra terms arising in an inflationary universe being multiplied by "slow-roll parameters" of inflationary dynamics. 

As first pointed out in [116] and later discussed in depth in [117, 118], it is crucial to measure the effects in terms of 
a physical clock as opposed to in terms of a background quantity which itself obtains corrections from the fluctuations. 
In the case of purely adiabatic fluctuations, it turns out that the effects of IR modes are not physically measurable 
[119, 120]. On the other hand, there arc effects of IR modes of entropy fluctuations which arc physically measurable 
[115]. 



Obviously, the expansion of space leads to a damping of these oscillations. 
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B. IR Effects 

In the previous subsection we have discussed infrared effects on cosmological fluctuations. What about infrared 
effects on the background? From the point of view of cosmological perturbations, in the same way that a mode with 
wavenumber k can combine with a mode with wavenumbcr k — k to cause a second order correction to the mode with 
wavenumber k, modes with wavenumbers k and —k lead to a second order correction to the cosmological background. 
We call this a "back-reaction" effect. 

In a long-standing research program, Tsamis and Woodard have been studying the back-reaction of infrared gravi- 
tational waves on the cosmological background [121]. They find a secular instability of the de Sitter background, on 
the basis of which they proposed a quantum gravitational model of inflation [122]. 

In [123], the back- reaction effect of scalar metric fluctuations in the dc Sitter phase of an inflationary cosmology 
was investigated (see [124] for a review of this method). The starting point was the following ansatz for metric and 
matter: 

g^(x,t) = g$(t) + g$(x,t) 
<p(x,t) = <pW(t) + <pW(x,t) (46) 

where g$(t) and <p^(t) form the cosmological background and g^J(x,t) and tp^(x,t) are the linear cosmological 
perturbations (the linearization parameter is the amplitude of the cosmological perturbations). 

This ansatz does not satisfy the Einstein equations beyond linear order. There are quadratic corrections to both 
the background and the fluctuations. Here we are interested in the corrections to the background (see [125] for an 
analysis of the back-reaction on the fluctuations using the same formalism) . We can absorb the quadratic corrections 

to the background metric into a new metric g^i (t) which equals the original background metric plus quadratic 
corrections which take into account that (46) does not solve the Einstein equations to quadratic order. 

To obtain the equations of motion for <?™g (i) we insert (46) into the full Einstein equations and expand to 
quadratic order. To extract the equation of motion for the modified background metric we take the spatial average of 
the resulting equations. There are terms quadratic in the linear metric fluctuations which we move to the right-hand 
side of the equations of motion to define an effective energy-momentum tensor for cosmological perturbations (see 

also [126, 127] where this method was first developed, albeit in a different context). The equation of motion for g^g 3 
which includes the quadratic corrections of first order fluctuations becomes 

= 8itG[T$+t^], (47) 

where G ^ v denotes the Einstein tensor, T^ v is the energy-momentum tensor of matter, and G stands for Newton's 
gravitational constant. 

In terms of the linear metric and matter fluctuations, the effective energy-momentum tensor takes the form 

^ = (^ L ^)' (48) 

where the first term contains the terms quadratic in the linear matter fluctuations, and the second term those quadratic 
in the linear metric perturbations. There are also products of metric and matter fluctuations which appear in the 
first term (sec [123] for the full expressions). 

The effective energy-momentum tensor obtains contributions from perturbations of all wavelengths. The con- 
tributions of all Fourier modes add up linearly. We can thus define an UV and an IR part which consist, respectively, 
of the contributions of sub-Hubble and super-Hubble modes. In the de Sitter limit of inflation, the UV part is time- 
independent (this follows both by symmetry and by explicit computation). On the other hand, due to the increasing 
phase space of IR modes, the IR part grows in time and hence dominates in the late time limit. 

We now sketch the evaluation of the IR part of Tn V . To be specific, we assume that matter takes the form of a 
scalar field ip. We work in longitudinal gauge in which the metric and matter including fluctuations take the form: 

ds 2 = a 2 [(l + 2$)d?7 2 - (1 - 2$)da; 2 ], tp = (p + Sip . (49) 

Note that the metric and matter fluctuations $ and Sip, respectively, are related by the Einstein constraint equations 



G nu 



(O.br) 



(50) 
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where V(<p) is the potential energy function of <p. To obtain the leading terms contributing to the IR part of tw, 
we can neglect spatial gradient terms and work to leading order in the slow-roll approximation for tp. With these 
approximations the IR part of t^ v takes the form of a cosmological constant with effective energy density 

/ h, =^ ( 2 ^-4vW>. (51) 



For simple inflationary models [e.g. V(tp) = ^m 2 ip 2 ] 

< 



and hence the back-reaction of IR modes takes the form of a negative contribution to the cosmological constant whose 
magnitude grows in time since the phase space of IR modes is increasing. 

We may conjecture that the back-reaction of IR modes thus leads to a dynamical relaxation mechanism for a bare 
positive cosmological constant A [124] (for a recent review on the cosmological constant problem, see Ref. [128]). We 
begin with a large bare A (in the presence of matter). The presence of A will lead to a period of inflation. This, in 
turn, will cause a sea of IR modes to build up. Then, back-reaction will set in. The effect of IR modes will lead to an 
effective cosmological constant which determines the full background metric which takes the form 

Aeff(t) = A-|p (br) W|, (52) 

and shows the onset of an instability. 

Naive extrapolation of the back-reaction effect to large times would lead to A e g(t) = at some time t = tiR. 
However, long before this happens the perturbative approach will break down. Assuming that the instability which 
is seen to leading order survives to a non-perturbative treatment, an interesting scenario emerges [124]: The back- 
reaction contribution ceases to increase in magnitude once Qa(£) < &m(t) (where fix is the contribution of the energy 
density of the substance X to the critical density) . since at that time the acceleration of space ceases and the infrared 
sea ceases to grow. The universe, however, is still expanding, and hence the matter energy density decreases. Thus, 
the effective cosmological constant once again raises its head. The upshot of this dynamics is that there is a dynamical 
fixed point with 

n A (t) ~ \ ■ (53) 

More precisely, oscillations of Q,\(t) about that value are expected. Thus, the back-reaction of long wavelength cosmo- 
logical perturbations has the potential of providing a dynamical relaxation mechanism for a large bare cosmological 
constant which explains dark energy and has no coincidence problem. 

Note that the above mechanism is not in conflict with causality since only modes are considered which are inside the 
horizon 8 . The effect can be locally described in terms of a time-dependent change in the spatial curvature constant 
k and A [129]. 

The key concern is that the above calculations have been performed as a function of a background time t. As 
conjectured in [116] and verified in [117, 118], in the case of purely adiabatic fluctuations (the energy densities of 
all components of matter being proportional) the leading IR back-reaction effect discussed above can be entirely 
absorbed by a second order time-reparametrization. In other words, if we calculate observables such as the local 
Hubble expansion rate as a function of a physical clock variable as opposed to the background time t, there is no 
leading order effect of the IR fluctuations. For adiabatic fluctuations, the only physical clock is the matter field <p, 
and we find [117] 

H 2 (tp,$) = H 2 (cp,0). (54) 

However, in the presence of multi-field systems with entropy fluctuations, it can be shown that the leading IR back- 
reaction effects are physically measurable [130] (see also [131]). If x is an entropy field (e.g. a field which represents 
the temperature of the CMB), then 

# 2 (X,$) + H 2 (x,0). (55) 

Thus, there appears to be evidence for an instability of the de Sitter phase of inflationary cosmology. What about de 
Sitter space itself? We will end with a very brief literature survey of work on this topic. 



Recall that in inflationary cosmology the horizon is larger than the Hubble radius by a factor which exponentially increases during 
inflation. 
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C. Stability of de Sitter 

De Sitter space-time is a classical solution of Einstein's field equations in the presence of a positive cosmological 
constant. According to the classical "no-hair" theorems [132, 133] the expanding branch is classically stable (the 
contracting branch is classically unstable because of the growth of fluctuations in such a phase). 

As discussed in the previous subsection, there is perturbative evidence that in the presence of matter and entropy 
fluctuations even the expanding branch of de Sitter space is semi-classically unstable. This question has been studied 
for a long time. Based on studies of the renormalization of the energy-momentum tensor of a massless scalar field 
in de Sitter space there have been claims of instability [134, 135]. There are also claims of instability due to particle 
production [136, 137], due to a thermodynamic instability [138] and due to a conformal anomaly [139]. The work of 
[121] mentioned in the previous subsection supports the claim of semi-classical instability. 

What about non-perturbative statements? In the context of four space-time dimensional gravity there have been 
long-standing claims that de Sitter space is unstable [140] (see also [141]). However, these claims have been disputed 
in [142]. 

There has been some interesting recent work on the stability issues of de Sitter space. In the context of three 
space-time dimensional pure gravity it has been shown that the method used to construct the partition function of 
quantum gravity which works in anti-de-Sittcr space breaks down in the presence of a positive cosmological constant 
[143] . In the context of the higher spin theory limit of string theory is has also been shown that de Sitter space does 
not seem to be a solution at the quantum level [144]. 

VI. CONCLUSIONS 

We have seen that the exponential expansion of space during the expanding branch of de Sitter space-time, and 
consequently also during an inflationary phase of an expanding FRWL universe, leads to important trans-Planckian 
issues. In particular, fluctuation modes which are being probed in current cosmological observations originate with 
wavelengths smaller than the Planck length at the onset of the phase of exponential expansion. Since the physics on 
trans-Planckian scales is unkown, this raises conceptual questions regarding to the robustness of the usual calculations 
in inflationary cosmology. 

One aspect of this problem, namely the fact that in an expanding space comoving modes continuously cross the 
boundary of the trans-Planckian zone of ignorance, is common to all expanding cosmologies. What is special to 
inflationary cosmology is that scales which are currently explored in cosmological observations emerge from this 
trans-Planckian sea. This is not the case in the two other cosmological paradigms mentioned earlier, namely the 
matter bounce scenario and string gas cosmology. 

These issues have been discussed in Section 2 of this article. Consequences for cosmological observations have 
been discussed in Section 3 and 4. Possibly the most important point is that due to the accelerated increase in the 
wavelength of scales, trans-Planckian physics can in fact be tested with current cosmological observations. We have 
seen that under the assumption that at some initial time all modes are in their local vacuum state, then oscillations 
in the spectrum of CMB anisotropics are induced. The strongest constraint on the amplitude of these oscillations 
comes from back-reaction considerations, namely from demanding that the produced particles do not destabilize the 
inflationary background. Such oscillations can be searched for in current CMB anisotropy data. With current results, 
a vanilla inflationary model without oscillations is still an excellent fit to the data. 

The accelerated expansion of space during an inflationary period leads to more conceptual problems. The increase 
of the wavelength of fixed comoving modes relative to the Hubble radius leads to a rapid increase of the phase space 
of such infrared (IR) modes. This may - in the presence of entropy fluctuations - lead to large IR effects, possibly 
even to a de-stabilization of the de Sitter phase. 

Acknowledgements 

We would like to thank C. Ringcval for having shared his numerical results with us. The work of R. B. is supported 
in part by an NSERC Discovery Grant and by funds from the Canada Research Chairs Program. J. M. would like to 
thank L. Sriramkumar for useful discussions and the University of Tokyo (Research Center for the Early Universe) 
where part of this work has been done. 



27 



References 



[1] A. Guth, "The Inflationary Universe: A Possible Solution To The Horizon And Flatness Problems," Phys. Rev. D 23, 
347 (1981). 

[2] R. Brout, F. Englert and E. Gunzig, "The Creation Of The Universe As A Quantum Phenomenon," Annals Phys. 115, 
78 (1978). 

[3] A. A. Starobinsky, "A New Type Of Isotropic Cosmological Models Without Singularity," Phys. Lett. B 91, 99 (1980). 
[4] K. Sato, "First Order Phase Transition Of A Vacuum And Expansion Of The Universe," Mon. Not. Roy. Astron. Soc. 
195, 467 (1981). 

[5] V. Mukhanov and G. Chibisov, "Quantum Fluctuation And Nonsingular Universe. (In Russian)," JETP Lett. 33, 532 

(1981) [Pisma Zh. Eksp. Teor. Fiz. 33, 549 (1981)]. 
[6] W. Press, "Spontaneous production of the Zel'dovich spectrum of cosmological fluctuations", Phys. Scr. 21, 702 (1980). 
[7] A. A. Starobinsky, "Spectrum of relict gravitational radiation and the early state of the universe," JETP Lett. 30, 682 

(1979) [Pisma Zh. Eksp. Teor. Fiz. 30, 719 (1979)]. 
[8] P. J. E. Peebles, J. T. Yu, "Primeval adiabatic perturbation in an expanding universe," Astrophys. J. 162, 815-836 (1970). 
[9] R. A. Sunyaev, Y. .B. Zeldovich, "Small scale fluctuations of relic radiation," Astrophys. Space Sci. 7, 3-19 (1970). 
[10] P. D. Mauskopf et al. [Boomerang Collaboration], "Measurement of a peak in the cosmic microwave background 

power spectrum from the North American test flight of BOOMERANG," Astrophys. J. 536, L59 (2000) [arXiv:astro- 

ph/9911444]. 

[11] C. L. Bennett et al., "First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and 

Basic Results," Astrophys. J. Suppl. 148, 1 (2003) [arXiv:astro-ph/0302207]. 
[12] R. H. Brandenberger and J. H. Kung, "Chaotic Inflation As An Attractor In Initial Condition Space," Phys. Rev. D 42, 

1008 (1990); 

R. H. Brandenberger, H. Feldman and J. Kung, "Initial conditions for chaotic inflation," Phys. Scripta T 36, 64 (1991). 
[13] V. F. Mukhanov, "Quantum Theory of Gauge Invariant Cosmological Perturbations," Sov. Phys. JETP 67, 1297 (1988) 

[Zh. Eksp. Teor. Fiz. 94N7, 1 (1988)]. 
[14] M. Sasaki, "Large Scale Quantum Fluctuations in the Inflationary Universe," Prog. Theor. Phys. 76, 1036 (1986). 
[15] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, "Theory of cosmological perturbations. Part 1. Classical 

perturbations. Part 2. Quantum theory of perturbations. Part 3. Extensions," Phys. Rept. 215, 203 (1992). 
[16] J. Martin, "Inflationary cosmological perturbations of quantum-mechanical origin", Lect. Notes Phys. 669, 199 (2005), 

[arXiv:hep-th/0406011] 

[17] R. H. Brandenberger, "Lectures on the theory of cosmological perturbations," Lect. Notes Phys. 646, 127 (2004) 
[arXiv:hep-th/0306071] . 

[18] R. H. Brandenberger, "Inflationary cosmology: Progress and problems," [arXiv:hep-ph/9910410]. 

[19] R. H. Brandenberger and J. Martin, "The Robustness of inflation to changes in superPlanck scale physics," 
Mod. Phys. Lett. A 16, 999 (2001), [arXiv:astro-ph/0005432]; 

J. Martin and R. H. Brandenberger, "The TransPlanckian problem of inflationary cosmology," Phys. Rev. D 63, 123501 
(2001), [arXiv:hep-th/0005209]. 

[20] A. Linde, V. Mukhanov and A. Vikman, "On adiabatic perturbations in the ekpyrotic scenario," JCAP 1002, 006 (2010) 
[arXiv:0912.0944 [hep-th]]. 

[21] J. C. Niemeyer, "Inflation with a Planck scale frequency cutoff," Phys. Rev. D 63, 123502 (2001) [astro-ph/0005533]. 
[22] W. G. Unruh, "Sonic analog of black holes and the effects of high frequencies on black hole evaporation," Phys. Rev. D 
51, 2827 (1995). 

[23] S. Corley and T. Jacobson, "Hawking spectrum and high frequency dispersion," Phys. Rev. D 54, 1568 (1996) [hep- 
th/9601073]. 

[24] P. Horava, "Quantum Gravity at a Lifshitz Point," Phys. Rev. D 79, 084008 (2009) [arXiv:0901.3775 [hep-th]]. 

[25] M. Lemoine, M. Lubo, J. Martin and J. -P. Uzan, "The Stress energy tensor for transPlanckian cosmology," Phys. Rev. 

D 65, 023510 (2002) [hep-th/0109128]. 
[26] R. H. Brandenberger and J. Martin, "On signatures of short distance physics in the cosmic microwave background," Int. 

J. Mod. Phys. A 17, 3663 (2002) [hep-th/0202142]. 
[27] R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, "A Generic estimate of transPlanckian modifications to the 

primordial power spectrum in inflation," Phys. Rev. D 66, 023518 (2002) [hep-th/0204129]; 

R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, "Imprints of short distance physics on inflationary cosmology," 
Phys. Rev. D 67, 063508 (2003) [hep-th/01 10226]; 

R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, "Inflation as a probe of short distance physics," Phys. Rev. D 64, 
103502 (2001) [hep-th/0104102]. 
[28] U. H. Danielsson, "A Note on inflation and transPlanckian physics," Phys. Rev. D 66, 023511 (2002) [hep-th/0203198]; 
U. H. Danielsson, "Inflation, holography, and the choice of vacuum in de Sitter space," JHEP 0207, 040 (2002) [hep- 
th/0205227]. 



28 



[29] J. C. Niemeyer, R. Parentani and D. Campo, "Minimal modifications of the primordial power spectrum from an adiabatic 

short distance cutoff," Phys. Rev. D 66, 083510 (2002) [hep-th/0206149]. 
[30] V. Bozza, M. Giovannini and G. Veneziano, "Cosmological perturbations from a new physics hypersurface," JCAP 0305, 

001 (2003) [hep-th/0302184]. 

[31] M. G. Jackson and K. Schalm, "Model-Independent Signatures of New Physics in non-Gaussianity," arXiv: 1202.0604 
[hep-th] ; 

M. G. Jackson and K. Schalm, "Model-Independent Signatures of New Physics in Slow-Roll Inflation," arXiv: 1104.0887 
[hep-th] ; 

M. G. Jackson and K. Schalm, "Model Independent Signatures of New Physics in the Inflationary Power Spectrum," 
Phys. Rev. Lett. 108, 111301 (2012) [arXiv: 1007.0185 [hep-th]]; 

K. Schalm, G. Shiu and J. P. van der Schaar, "The Cosmological vacuum ambiguity, effective actions, and transplanckian 
effects in inflation," AIP Conf. Proc. 743, 362 (2005) [hep-th/0412288]; 

K. Schalm, G. Shiu and J. P. van der Schaar, "Decoupling in an expanding universe: Boundary RG flow affects initial 

conditions for inflation," JHEP 0404, 076 (2004) [hep-th/0401164]. 
[32] J. Martin and R. Brandenberger, "On the dependence of the spectra of fluctuations in inflationary cosmology on trans- 

Planckian physics," Phys. Rev. D 68, 063513 (2003) [hep-th/0305161]. 
[33] T. S. Bunch and P. C. W. Davies, "Quantum Field Theory in de Sitter Space: Renormalization by Point Splitting," Proc. 

Roy. Soc. Lond. A 360, 117 (1978). 
[34] E. A. Tagirov, "Consequences of field quantization in de Sitter type cosmological models," Annals Phys. 76, 561 (1973); 

N. A. Chernikov and E. A. Tagirov, "Quantum theory of scalar fields in de Sitter space-time," Annales Poincare Phys. 

Theor. A 9, 109 (1968); 

E. Mottola, "Particle Creation in de Sitter Space," Phys. Rev. D 31, 754 (1985); 

B. Allen, "Vacuum States in de Sitter Space," Phys. Rev. D 32, 3136 (1985). 

[35] K. Goldstein and D. A. Lowe, "Initial state effects on the cosmic microwave background and transPlanckian physics," 

Phys. Rev. D 67, 063502 (2003) [hep-th/0208167]. 
[36] N. Kaloper, M. Kleban, A. Lawrence, S. Shenker and L. Susskind, "Initial conditions for inflation," JHEP 0211, 037 

(2002) [hep-th/0209231]. 

[37] M. B. Einhorn and F. Larsen, "Interacting quantum field theory in de Sitter vacua," Phys. Rev. D 67, 024001 (2003) 
[hep-th/0209159]. 

[38] K. Goldstein and D. A. Lowe, "A Note on alpha vacua and interacting field theory in de Sitter space," Nucl. Phys. B 

669, 325 (2003) [hep-th/0302050] . 
[39] A. Borde and A. Vilenkin, "Eternal inflation and the initial singularity," Phys. Rev. Lett. 72, 3305 (1994) [gr-qc/9312022]. 
[40] N. Kaloper, M. Kleban, A. E. Lawrence and S. Shenker, "Signatures of short distance physics in the cosmic microwave 

background," Phys. Rev. D 66, 123510 (2002) [hep-th/0201158]. 
[41] C. P. Burgess, J. M. Cline, F. Lemieux and R. Holman, "Are inflationary predictions sensitive to very high-energy 

physics?," JHEP 0302, 048 (2003) [hep-th/0210233]; 

C. P. Burgess, J. M. Cline and R. Holman, "Effective field theories and inflation," JCAP 0310, 004 (2003) [hep- 
th/0306079]. 

[42] J. Martin, "Inflationary perturbations: The Cosmological Schwinger effect", Lect. Notes Phys. 738, 193 (2008), 
[arXiv:0704.3540]. 

[43] S. F. Hassan and M. S. Sloth, "TransPlanckian effects in inflationary cosmology and the modified uncertainty principle," 

Nucl. Phys. B 674, 434 (2003) [hep-th/0204110]. 
[44] C. -S. Chu, B. R. Greene and G. Shiu, "Remarks on inflation and noncommutative geometry," Mod. Phys. Lett. A 16, 

2231 (2001) [hep-th/0011241]. 

[45] R. Brandenberger and P. -M. Ho, "Noncommutative space-time, stringy space-time uncertainty principle, and density 

fluctuations," Phys. Rev. D 66, 023517 (2002) [AAPPS Bull. 12N1, 10 (2002)] [hep-th/0203119]. 
[46] L. Senatore, private communication. 

[47] Y. F. Cai, T. Qiu, R. Brandenberger, Y. S. Piao and X. Zhang, "On Perturbations of Quintom Bounce," JCAP 0803, 

013 (2008) [arXiv:0711.2187 [hep-th]]. 
[48] I. Florakis, C. Kounnas, H. Partouche and N. Toumbas, "Non-singular string cosmology in a 2d Hybrid model," Nucl. 

Phys. B 844, 89 (2011) [arXiv:1008.5129 [hep-th]]; 

C. Kounnas, H. Partouche and N. Toumbas, "Thermal duality and non-singular cosmology in d-dimensional superstrings," 
Nucl. Phys. B 855, 280 (2012) [arXiv: 1106.0946 [hep-th]]; 

C. Kounnas, H. Partouche and N. Toumbas, "S-brane to thermal non-singular string cosmology," Class. Quant. Grav. 
29, 095014 (2012) [arXiv:1111.5816 [hep-th]]. 
[49] R. Brandenberger, "Matter Bounce in Horava-Lifshitz Cosmology," Phys. Rev. D 80, 043516 (2009) [arXiv:0904.2835 
[hep-th]]. 

[50] R. Brandenberger and X. -m. Zhang, "The Trans-Planckian Problem for Inflationary Cosmology Revisited," 
arXiv:0903.2065 [hep-th]. 

[51] D. Wands, "Duality invariance of cosmological perturbation spectra," Phys. Rev. D 60, 023507 (1999) [arXiv:gr- 
qc/9809062]. 

[52] F. Finelli and R. Brandenberger, "On the generation of a scale-invariant spectrum of adiabatic fluctuations in cosmological 

models with a contracting phase," Phys. Rev. D 65, 103522 (2002) [arXiv:hep-th/0112249]. 
[53] S. Mukohyama, "Horava-Lifshitz Cosmology: A Review," Class. Quant. Grav. 27, 223101 (2010) [arXiv:1007.5199 [hep- 



29 



th]]. 

[54] A. Cerioni and R. H. Brandenberger, "Cosmological Perturbations in the Projectable Version of Horava-Lifshitz Gravity," 

JCAP 1108, 015 (2011) [arXiv: 1007. 1006 [hep-th]]. 
[55] A. Cerioni and R. H. Brandenberger, "Cosmological Perturbations in the 'Healthy Extension' of Horava-Lifshitz gravity," 

arXiv: 1008.3589 [hep-th]. 

[56] E. G. M. Ferreira and R. Brandenberger, "The Trans-Planckian Problem in the Healthy Extension of Horava-Lifshitz 

Gravity," Phys. Rev. D 86, 043514 (2012) [arXiv: 1204.5239 [hep-th]]. 
[57] R. H. Brandenberger and C. Vafa, "Superstrings in the Early Universe," Nucl. Phys. B 316, 391 (1989). 
[58] J. Kripfganz and H. Perlt, "Cosmological Impact Of Winding Strings," Class. Quant. Grav. 5, 453 (1988). 
[59] A. Nayeri, R. H. Brandenberger and C. Vafa, "Producing a scale-invariant spectrum of perturbations in a Hagedorn phase 

of string cosmology," Phys. Rev. Lett. 97, 021302 (2006) [arXiv:hep-th/0511140]. 
[60] R. H. Brandenberger, A. Nayeri, S. P. Patil and C. Vafa, "String gas cosmology and structure formation," Int. J. Mod. 

Phys. A 22, 3621 (2007) [arXiv:hep-th/0608121]. 
[61] R. H. Brandenberger, "Challenges for string gas cosmology," arXiv:hep-th/0509099; 

T. Battefeld and S. Watson, "String gas cosmology," Rev. Mod. Phys. 78, 435 (2006) [arXiv:hep-th/0510022]. 
[62] R. Hagedorn, "Statistical Thermodynamics Of Strong Interactions At High- Energies," Nuovo Cim. Suppl. 3, 147 (1965). 
[63] R. H. Brandenberger, A. Nayeri, S. P. Patil and C. Vafa, "Tensor modes from a primordial Hagedorn phase of string 

cosmology," Phys. Rev. Lett. 98, 231302 (2007) [arXiv:hep-th/0604126]. 
[64] N. Weiss, "Constraints On Hamiltonian Lattice Formulations Of Field Theories In An Expanding Universe," Phys. Rev. 

D 32, 3228 (1985). 

[65] R. H. Brandenberger, "The Matter Bounce Alternative to Inflationary Cosmology," arXiv: 1206.4196 [astro-ph.CO]; 

R. H. Brandenberger, "Introduction to Early Universe Cosmology," PoS ICFI 2010, 001 (2010) [arXiv: 1103.2271 [astro- 
ph.CO]]. 

[66] H. Li, J. Q. Xia, R. Brandenberger and X. Zhang, "Constraints on Models with a Break in the Primordial Power 

Spectrum," arXiv:0903.3725 [astro-ph.CO]. 
[67] Y. F. Cai, W. Xue, R. Brandenberger and X. Zhang, "Non-Gaussianity in a Matter Bounce," JCAP 0905, 011 (2009) 

[arXiv:0903.0631 [astro-ph.CO]]. 
[68] X. Chen, "Primordial Features as Evidence for Inflation", JCAP 1201, 038 (2012), [arXiv: 1104. 1323]. 
[69] X. Chen, "Searching for Standard Clocks in the Primordial Universe", JCAP 1208, 014 (2012), [arXiv:1205.6085]. 
[70] A. A. Starobinsky, "Robustness of the inflationary perturbation spectrum to transPlanckian physics," Pisma Zh. Eksp. 

Teor. Fiz. 73, 415 (2001) [JETP Lett. 73, 371 (2001)] [astro-ph/0104043]; 

A. A. Starobinsky and I. I. Tkachev, "Trans-Planckian particle creation in cosmology and ultra-high energy cosmic rays," 
JETP Lett. 76, 235 (2002) [Pisma Zh. Eksp. Teor. Fiz. 76, 291 (2002)] [astro- ph/02075 72]. 
[71] T. Tanaka, "A Comment on transPlanckian physics in inflationary universe," astro-ph/00 12431. 

[72] R. H. Brandenberger and J. Martin, Back-reaction and the trans-Planckian problem of inflation revisited", 

Phys. Rev. D 71, 023504 (2005), [arXiv:hep-th/0410223]. 
[73] U. Danielsson, " Transplanckian energy production and slow roll inflation", Phys. Rev. D 71, 023516 (2005), [arXiv:hep- 

th/0411172]. 

[74] J. Martin and C. Ringeval, "Superimposed oscillations in the WMAP data?", Phys. Rev. D 69, 083515 (2004), 
[arXiv:astro-ph/0310382]. 

[75] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, M. Nolta et al, "Seven- Year Wilkinson Microwave Anisotropy 
Probe (WMAP) Observations: Power Spectra and WMAP-Derived Parameters", Astrophys. J. Suppl. 192, 16 (2011), 
[arXiv:1001.4635]. 

[76] E. Komatsu et al, "Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpre- 
tation", Astrophys. J. Suppl. 192, 18 (2011), [arXiv:1001.4538]. 

[77] J. Martin and C. Ringeval, "Exploring the superimposed oscillations parameter space", JCAP 0501, 007 (2005), 
[arXiv:hep-ph/0405249] . 

[78] J. Martin and C. Ringeval, "Addendum to 'Superimposed oscillations in the WMAP data?", Phys. Rev. D 69, 127303 

(2004) , [arXiv:astro-ph/0402609] . 
[79] J. Martin and C. Ringeval, "Inflation after WMAP3: Confronting the Slow- Roll and Exact Power Spectra to CMB Data", 

JCAP 0608, 009 (2006), [arXiv:astro-ph/0605367]. 
[80] L. Lorenz, J. Martin and C. Ringeval, "Brane inflation and the WMAP data: A Bayesian analysis", JCAP 0804, 001 

(2008), [arXiv:0709.3758]. 

[81] L. Lorenz, J. Martin and C. Ringeval, "Constraints on Kinetically Modified Inflation from WMAP5", Phys. Rev. D 78, 

063543 (2008), [arXiv:0807.2414]. 
[82] J. Martin and C. Ringeval, "First CMB Constraints on the Inflationary Reheating Temperature", Phys. Rev. D 82, 

023511 (2010), [arXiv: 1004.5525]. 
[83] J. Martin, C. Ringeval and R. Trotta, "Hunting Down the Best Model of Inflation with Bayesian Evidence", 

Phys. Rev. D 83, 063524 (2011), [arXiv:1009.4157]. 
[84] A. Lewis, A. Challinor and A. Lasenby, "Efficient computation of CMB anisotropies in closed FRW models" , Astrophys. J. 

538, 473 (2000), [arXiv:astro-ph/astro-ph/9911177]. 
[85] A. Lewis and S. Bridle, "Cosmological parameters from CMB and other data: A Monte Carlo approach", Phys. Rev. D 

66, 103511 (2002), [arXiv:astro-ph/astro-ph/0205436]. 
[86] C. Ringeval, private communication. 



30 



[87] R. Easther, W. Kinney and H. Peiris, "Observing trans-Planckian signatures in the cosmic microwave background", JCAP 

0505, 009 (2005), [arXiv:astro-ph/0412613]. 
[88] A. Gangui, F. Lucchin, S. Matarrese and S. Mollerach, "The Three point correlation function of the cosmic microwave 

background in inflationary models", Astrophys. J. 430, 447 (1994), [arXiv:astro-ph/9312033]. 
[89] A. Gangui, "NonGaussian effects in the cosmic microwave background from inflation", Phys. Rev. D 50, 3684 (1994), 

[arXiv:astro-ph/9406014]. 

[90] L. Wang and M. Kamionkowski, "The Cosmic microwave background bispectrum and inflation", Phys. Rev. D 61, 063504 

(2000) , [arXiv:astro-ph/9907431] . 
[91] A. Gangui and J. Martin, "Cosmic microwave background bispectrum and slow roll inflation", Mon. Not. Roy. Astron. Soc. 

313, 323 (1999)", [arXiv:astro-ph/9908009]. 
[92] A. Gangui and J. Martin, "Best unbiased estimators for the three point correlators of the cosmic microwave background 

radiation", Phys. Rev. D 62, 103004 (2000), [arXiv:astro-ph/0001361]. 
[93] J. M. Maldacena, "Non-Gaussian features of primordial fluctuations in single field inflationary models", JHEP 0305, 013 

(2003), [arXiv:astro-ph/0210603]. 
[94] J. Martin, A. Riazuelo and M. Sakellariadou, "Nonvacuum initial states for cosmological perturbations of quantum 

mechanical origin", Phys. Rev. D 61, 083518 (2000), [arXiv:astro-ph/9904167]. 
[95] A. Gangui, J. Martin and M. Sakellariadou, "Single field inflation and non-Gaussianity" , Phys. Rev. D 66, 083502 (2002), 

[arXiv:astro-ph/0205202] . 

[96] A. Gangui, J. Martin and M. Sakellariadou, "Non- vacuum initial states for the cosmological perturbations and the back 

reaction problem of inflation", [arXiv:astro-ph/0206023]. 
[97] R. Holman and A. Tolley, "Enhanced Non-Gaussianity from Excited Initial States", JCAP 0805, 001 (2008), 

[arXiv:0710.1302]. 

[98] P. Meerburg, J. P. van der Schaar and P. S. Corasaniti, " Signatures of Initial State Modifications on Bispectrum Statistics" , 

JCAP 0905, 018 (2009), [arXiv:0901.4044]. 
[99] I. Agullo and L. Parker, " Non-gaussianities and the Stimulated creation of quanta in the inflationary universe", 

Phys. Rev. D 83, 063526 (2011), [arXiv:1010.5766]. 
[100] J. Ganc, "Calculating the local-type fNL for slow-roll inflation with a non- vacuum initial state", Phys. Rev. D 84, 063514 

(2011), [arXiv:1104.0244]. 

[101] D. Chialva, "Signatures of very high energy physics in the squeezed limit of the bispectrum (violation of Maldacena's 
condition)", [arXiv:1108.4203]. 

[102] X. Chen, M. Huang, K. Shamit and G. Shiu, "Observational signatures and non-Gaussianities of general single field 

inflation", JCAP 0701, 002 (2007), [arXiv:hep-th/0605045] . 
[103] X. Cheng, "Primordial Non-Gaussianities from Inflation Models", Adv. Astron. 2010, 638979 (2010), [arXiv:1002.1416]. 
[104] D. Seery and J. E. Lidsey, "Primordial non-Gaussianities in single field inflation", JCAP 0506, 003 (2005), [arXiv:astro- 

ph/0503692]. 

[105] X. Cheng, R. Easther and E. Lim, "Generation and Characterization of Large Non-Gaussianities in Single Field Inflation", 

JCAP 0804, 010 (2008), [arXiv:0801.3295]. 
[106] X. Cheng, R. Easther and E. Lim, "Large Non-Gaussianities in Single Field Inflation", JCAP 0706, 023 (2007), 

[arXiv:astro-ph/0611645]. 

[107] J. Martin and L. Sriramkumar, "The scalar bi-spectrum in the Starobinsky model: The equilateral case", JCAP 1201, 

008 (2012), [arXiv:1109.5838]. 
[108] D. Hazra, L. Sriramkumar and J. Martin, "On the discriminating power of /nl", [arXiv:1201.0926]. 

[109] D. Hazra, J. Martin and L. Sriramkumar, "The scalar bi-spectrum during preheating in single field inflationary models", 

Phys. Rev. D 86, 063523 (2012), [arXiv:1206.0442]. 
[110] S. Renaux-Petel, "On the squeezed limit of the bispectrum in general single field inflation", JCAP 1010, 020 (2010), 

[arXiv: 1008.0260]. 

[Ill] J. Martin, H. Motohashi and T. Suyama, "Ultra Slow-Roll Inflation and the non-Gaussianity Consistency Relation", 
[arXiv:1211.0083]. 

[112] J. Ganc and E. Komatsu, "Scale-dependent bias of galaxies and mu-type distortion of the cosmic microwave background 
spectrum from single-field inflation with a modified initial state", Phys. Rev. D 86, 023518 (2012), [arXiv:1204.4241]. 

[113] W. Xue, K. Dasgupta and R. Brandenberger, "Cosmological UV/IR Divergences and de-Sitter Spacetime," Phys. Rev. 
D 83, 083520 (2011) [arXiv:1103.0285 [hep-th]]. 

[114] S. Weinberg, "Quantum contributions to cosmological correlations," Phys. Rev. D 72, 043514 (2005) [arXiv:hep- 
th/0506236]; 

S. B. Giddings, M. S. Sloth, "Semiclassical relations and IR effects in de Sitter and slow-roll space-times," JCAP 1101, 
023 (2011). [arXiv:1005.1056 [hep-th]]; 

S. B. Giddings and M. S. Sloth, "Cosmological observables, IR growth of fluctuations, and scale-dependent anisotropies," 
Phys. Rev. D 84, 063528 (2011) [arXiv:1104.0002 [hep-th]]; 

D. Seery, "Infrared effects in inflationary correlation functions," Class. Quant. Grav. 27, 124005 (2010). [arXiv: 1005. 1649 
[astro-ph.CO]]; 

M. S. Sloth, "On the one loop corrections to inflation and the CMB anisotropies," Nucl. Phys. B748, 149-169 (2006). 
[astro-ph/0604488]; 

M. S. Sloth, "On the one loop corrections to inflation. II. The Consistency relation," Nucl. Phys. B775, 78-94 (2007). 
[hep-th/0612138]; 



31 



D. Seery, "One-loop corrections to a scalar field during inflation," JCAP 0711, 025 (2007) [arXiv:0707.3377 [astro-ph]]; 

D. Seery, "One-loop corrections to the curvature perturbation from inflation," JCAP 0802, 006 (2008). [arXiv:0707.3378 
[astro-ph]]. N. Bartolo, A. Riotto, "Possibly Large Corrections to the Inflationary Observables," Mod. Phys. Lett. A23, 
857-862 (2008) [arXiv:0711.4003 [astro-ph]]; 

N. Bartolo, S. Matarrese, M. Pietroni et al, "On the Physical Significance of Infra-red Corrections to Inflationary Ob- 
servables," JCAP 0801, 015 (2008) [arXiv:0711.4263 [astro-ph]]; 

Y. Urakawa, K. -i. Maeda, "One-loop Corrections to Scalar and Tensor Perturbations during Inflation in Stochastic Grav- 
ity," Phys. Rev. D78, 064004 (2008). [arXiv:0801.0126 [hep-th]]; 

K. Enqvist, S. Nurmi, D. Podolsky et al, "On the divergences of inflationary superhorizon perturbations," JCAP 0804, 
025 (2008). [arXiv:0802.0395 [astro-ph]]; 

E. Dimastrogiovanni, N. Bartolo, "One-loop graviton corrections to the curvature perturbation from inflation," JCAP 
0811, 016 (2008); [arXiv:0807.2790 [astro-ph]]; 

X. Gao, F. Xu, "Loop Corrections to Cosmological Perturbations in Multi-field Inflationary Models," JCAP 0907, 042 

(2009) . [arXiv:0905.0405 [hep-th]]; 

D. Campo, "Quantum corrections during inflation and conservation of adiabatic perturbations," Phys. Rev. D81, 043535 

(2010) . [arXiv:0908.3642 [hep-th]]; 

J. Kumar, L. Leblond, A. Rajaraman, "Scale Dependent Local Non-Gaussianity from Loops," JCAP 1004, 024 (2010). 
[arXiv:0909.2040 [astro-ph. CO]]. N. Bartolo, E. Dimastrogiovanni, A. Vallinotto, "One-loop corrections to the power spec- 
trum in general single-field inflation," JCAP 1011, 003 (2010). [arXiv:1006.0196 [astro-ph. CO]]; 

C. T. Byrnes, M. Gerstenlauer, A. Hebecker, S. Nurmi and G. Tasinato, "Inflationary Infrared Divergences: Geometry 
of the Reheating Surface vs. delta N Formalism," JCAP 1008, 006 (2010) [arXiv: 1005.3307 [hep-th]]; 

E. O. Kahya, V. K. Onemli and R. P. Woodard, "The Zeta-Zeta Correlator Is Time Dependent," Phys. Lett. B 694, 101 
(2010) [arXiv:1006.3999 [astro-ph.CO]]; 

S. P. Miao and R. P. Woodard, "Issues Concerning Loop Corrections to the Primordial Power Spectra," JCAP 1207, 008 
(2012) [arXiv:1204.1784 [astro-ph.CO]]; 

S. P. Miao, N. C. Tsamis and R. P. Woodard, "De Sitter Breaking through Infrared Divergences," J. Math. Phys. 51, 

072503 (2010) [arXiv:1002.4037 [gr-qc]]. 
[115] W. Xue, X. Gao and R. Brandenberger, "IR Divergences in Inflation and Entropy Perturbations," JCAP 1206, 035 

(2012) [arXiv: 120 1.0768 [hep-th]]. 
[116] W. Unruh, "Cosmological long wavelength perturbations," astro-ph/9802323. 

[117] G. Geshnizjani and R. Brandenberger, "Back reaction and local cosmological expansion rate," Phys. Rev. D 66, 123507 
(2002) [gr-qc/0204074]. 

[118] L. R. Abramo and R. P. Woodard, "No one loop back-reaction in chaotic inflation," Phys. Rev. D 65, 063515 (2002) 
[arXiv:astro-ph/0109272]. 

[119] M. Gerstenlauer, A. Hebecker and G. Tasinato, "Inflationary Correlation Functions without Infrared Divergences," 
arXiv:1102.0560 [astro-ph.CO]; 

Y. Urakawa, T. Tanaka, "IR divergence does not affect the gauge-invariant curvature perturbation," Phys. Rev. D82, 
121301 (2010). [arXiv: 1007.0468 [hep-th]]; 

Y. Urakawa, T. Tanaka, "Influence on Observation from IR Divergence during Inflation. I.," Prog. Theor. Phys. 122, 

779-803 (2009). [arXiv:0902.3209 [hep-th]]. 
[120] L. Senatore and M. Zaldarriaga, "On Loops in Inflation," JHEP 1012, 008 (2010) [arXiv:0912.2734 [hep-th]]; 

L. Senatore and M. Zaldarriaga, "On Loops in Inflation II: IR Effects in Single Clock Inflation," arXiv: 1203. 6354 [hep-th]; 

G. L. Pimentel, L. Senatore and M. Zaldarriaga, "On Loops in Inflation III: Time Independence of zeta in Single Clock 

Inflation," JHEP 1207, 166 (2012) [arXiv: 1203.6651 [hep-th]]. 
[121] N. C. Tsamis and R. P. Woodard, "Relaxing The Cosmological Constant," Phys. Lett. B 301, 351 (1993); 

N. C. Tsamis and R. P. Woodard, "The Physical basis for infrared divergences in inflationary quantum gravity," Class. 

Quant. Grav. 11, 2969 (1994); 

N. C. Tsamis and R. P. Woodard, "Strong infrared effects in quantum gravity," Annals Phys. 238, 1 (1995); 

N. C. Tsamis and R. P. Woodard, "The quantum gravitational back-reaction on inflation," Annals Phys. 253, 1 (1997) 

[arXiv:hep-ph/9602316]; 

N. C. Tsamis and R. P. Woodard, "Quantum Gravity Slows Inflation," Nucl. Phys. B 474, 235 (1996) [arXiv:hep- 
ph/9602315]; 

N. C. Tsamis and R. P. Woodard, "One Loop Graviton Self-Energy In A Locally De Sitter Background," Phys. Rev. D 
54, 2621 (1996) [arXiv:hep-ph/9602317]. 
[122] N. C. Tsamis and R. P. Woodard, "A Phenomenological Model for the Early Universe," Phys. Rev. D 80, 083512 (2009) 
[arXiv:0904.2368 [gr-qc]]. 

[123] L. R. W. Abramo, R. H. Brandenberger and V. F. Mukhanov, "The energy-momentum tensor for cosmological pertur- 
bations," Phys. Rev. D 56, 3248 (1997) [arXiv:gr-qc/9704037]; 

V. F. Mukhanov, L. R. W. Abramo and R. H. Brandenberger, "On the back reaction problem for gravitational pertur- 
bations," Phys. Rev. Lett. 78, 1624 (1997) [arXiv:gr-qc/9609026] ; 

L. R. W. Abramo and R. P. Woodard, "One loop back reaction on chaotic inflation," Phys. Rev. D 60, 044010 (1999) 
[arXiv: astro-ph/98 1 1430] . 

[124] R. H. Brandenberger, "Back reaction of cosmological perturbations and the cosmological constant problem," hep- 
th/0210165. 



32 



[125] P. Martineau and R. H. Brandenberger, "The Effects of gravitational back-reaction on cosmological perturbations," Phys. 
Rev. D 72, 023507 (2005) [astro- ph/0505236]. 

[126] R. A. Isaacson, "Gravitational Radiation in the Limit of High Frequency. II. Nonlinear Terms and the Ef fective Stress 
Tensor," Phys. Rev. 166, 1272 (1968). 

[127] D. R. Brill and J. B. Hartle, "Method of the Self-Consistent Field in General Relativity and its Application to the 
Gravitational Geon," Phys. Rev. 135, B271 (1964). 

[128] J. Martin, "Everything You Always Wanted To Know About The Cosmological Constant Problem (But Were Afraid To 
Ask)", Comptes Rendus Physique 13, 566 (2012), [arXiv: 1205.3365]. 

[129] R. H. Brandenberger and C. S. Lam, "Back-reaction of cosmological perturbations in the infinite wavelength approxima- 
tion," hep-th/0407048. 

[130] G. Geshnizjani and R. Brandenberger, "Back reaction of perturbations in two scalar field inflationary models," JCAP 

0504, 006 (2005) [hep-th/0310265]. 
[131] Y. Urakawa, "Influence of gauge artifact on adiabatic and entropy perturbations during inflation," Prog. Theor. Phys. 

126, 961 (2011) [arXiv:1105.1078 [hep-th]]. 
[132] A. A. Starobinsky, "Isotropization of arbitrary cosmological expansion given an effective cosmological constant," JETP 

Lett. 37, 66 (1983). 

[133] R. M. Wald, "Asymptotic behavior of homogeneous cosmological models in the presence of a positive cosmological 

constant," Phys. Rev. D 28, 2118 (1983). 
[134] L. H. Ford, "Quantum Instability Of De Sitter Space-Time," Phys. Rev. D 31, 710 (1985). 

[135] J. H. Traschen and C. T. Hill, "Instability Of De Sitter Space On Short Time Scales," Phys. Rev. D 33, 3519 (1986). 
[136] N. P. Myhrvold, "Runaway Particle Production In De Sitter Space," Phys. Rev. D 28, 2439 (1983). 
[137] E. Mottola, "Particle Creation In De Sitter Space," Phys. Rev. D 31, 754 (1985); 

P. Mazur and E. Mottola, "Spontaneous Breaking Of De Sitter Symmetry By Radiative Effects," Nucl. Phys. B 278, 694 

(1986). 

[138] E. Mottola, "Thermodynamic Instability Of De Sitter Space," Phys. Rev. D 33, 1616 (1986). 

[139] I. Antoniadis and E. Mottola, "4-D quantum gravity in the conformal sector," Phys. Rev. D 45, 2013 (1992). 

[140] A. M. Polyakov, "De Sitter Space and Eternity," Nucl. Phys. B 797, 199 (2008) [arXiv:0709.2899 [hep-th]]; 

E. T. Akhmedov, P. V. Buividovich and D. A. Singleton, "De Sitter space and perpetuum mobile," arXiv:0905.2742 

[gr-qc]; 

A. M. Polyakov, "Decay of Vacuum Energy," Nucl. Phys. B 834, 316 (2010) [arXiv:0912.5503 [hep-th]]; 

D. Krotov and A. M. Polyakov, "Infrared Sensitivity of Unstable Vacua," arXiv: 1012.2107 [hep-th]; 
A. M. Polyakov, "Infrared instability of the de Sitter space," arXiv:1209.4135 [hep-th]. 

[141] E. T. Akhmedov, "IR divergences and kinetic equation in de Sitter space. Poincare patch: Principal series," 
arXiv:1110.2257 [hep-th]; 

K. Larjo and D. A. Lowe, "Initial states and infrared physics in locally de Sitter spacetime," arXiv: 11 12.5425 [hep-th]; 

E. T. Akhmedov, "On the physical meaning and consequences of the loop IR divergences in global dS," arXiv: 1209.4448 
[hep-th]. 

[142] A. Rajaraman, J. Kumar, L. Leblond, "Constructing Infrared Finite Propagators in Inflating Space-time," Phys. Rev. 
D82, 023525 (2010). [arXiv: 1002.4214 [hep-th]]; 

D. Marolf and I. A. Morrison, "The IR stability of de Sitter: Loop corrections to scalar propagators," Phys. Rev. D 82, 
105032 (2010) [arXiv:1006.0035 [gr-qc]]; 

A. Rajaraman, "On the proper treatment of massless fields in Euclidean de Sitter space," Phys. Rev. D82, 123522 (2010). 
[arXiv: 1008. 1271 [hep-th]]; 

D. Marolf and I. A. Morrison, "The IR stability of de Sitter QFT: results at all orders," arXiv: 1010.5327 [gr-qc]; 

A. Higuchi, D. Marolf and I. A. Morrison, "On the Equivalence between Euclidean and In-In Formalisms in de Sitter 

QFT," arXiv: 1012.3415 [gr-qc]. 

[143] A. Castro, N. Lashkari and A. Maloney, "A de Sitter Farey Tail," Phys. Rev. D 83, 124027 (2011) [arXiv: 1103.4620 
[hep-th]]; 

A. Castro and A. Maloney, "The Wave Function of Quantum de Sitter," arXiv: 1209. 5757 [hep-th]. 
[144] D. Anninos, F. Denef and D. Harlow, "The Wave Function of Vasiliev's Universe - A Few Slices Thereof," arXiv: 1207.5517 
[hep-th]. 



